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consideration was given to ground state molecular structures, of low molecular weight, which exhibit a 
high positive heat of formation. Calculations to date have been carried out on: 1) light element C} v 
and C 2 v structures (H 4 , U 3 H, UH 3 , U 4 ). 2) azide-like structures (FN' 3 , a-N 202 . FNCO CO 3 . HN 3 , 
FNBF); 3) cyclic boron structures (B 3 H 3 , B 2 H 2 NH); and 4) hypervalent structures (NF 5 , PF 5 ). 

We find that the ground states of H 4 (tetrahydrogen), UH 3 and U 4 . as C 3 V structures, are unsta¬ 
ble toward dissociation into diatomic fragments. These molecules do not look promising as storable, 
high energy chemical species. We have studied the U 3 H system in both C 3 V (trigonal pyramid) and 
C 2 v (planar kite-like) symmetries. Although both structures are vibrationally stable, with the C 2 v 
symmetry slightly lower in energy at the MP2/6-311G** level of theory, an analysis of the bond 
character in the C 3 V and C 2 v structures indicate that only the latter is thermodynamically stable. Using 
either O 2 or F 2 as an oxidizer, U 3 H appears to give about 10% improvement in I sp over LiH. The 
optimum I sp for the U 3 H/H 2 /O 2 system is 456 seconds, which can be compared with 441 seconds for 
the LiH/H 2 /C >2 system and 441 seconds for the H 2 /O 2 system run at P c = 800 psia, A/A* = 20:1. 

All of the azide-like structures are vibrationally stable with high (400-500 kJ/mol) positive heats 
of formation. The best candidate appears to be asymmetric dinitrogen dioxide. a-N 2 C> 2 . which exhib 
its a heat of formation of +500 kJ/mol and can deliver 450 seconds I sp as a hybrid propellant with 
hydrogen, and advanced performance in solid propellant formulations as a replacement for AP The 
energetics and pathways for decomposition of FN r 3 have also been studied. Our calculations, carried 
out to compliment the experimental studies of FN 3 at Rockwell Science Center, indicate that 
NFJa 1 A] should be efficiently produced by the thermal decomposition of FN 3 . This shows promise as 
a method to produce a short wavelength chemical laser. The structure of this molecule, its predicted 
vibrational spectrum and the role of low-lying triplet states have been examined. 

None of the cyclic boron structures examined to date, which include B 3 H 3 and B 2 H 2 NH, are 
highly energetic. Their chemistry is further complicated by condensation reactions which cause con¬ 
version to higher molecular weight structures. 

The hypervalent NF 5 molecule has been studied as a D 3 h structure, as the ionic NF 4 *F~ form and 
as an inverted C 3 V structure. Detailed calculations indicate stability for NF 5 only as a D 3 h structure 
Correlated energy calculations were carried out for NF 5 and for several azide structures Our analysis 
indicates that NF 5 is a thermodynamically stable compound with a predicted heat of formation AHf° 
(0 K) = + 25.1 kJ/mol. An H 2 /F 2 /NF 5 propellant system is predicted to have an optimum I sp of 466 
seconds (P c = 800 psia, A/A* = 20:1). yielding no significant improvement over the H 2 /F 2 system. 
However, if the physical properties of NF 5 indicate the possibility of non-cryogenic storage, NF 5 may 
be a promising candidate for an earth storable oxidizer. 
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SUMMARY 


A theoretical research investigation of light element species that may have potential for use 
as advanced propellants or as energy sources in a chemical propulsion system has been undertak¬ 
en. Included in this research program were detailed quantum mechanical studies of the electronic 
struct'_ ,r e of such species, an analysis of their chemical (and thermodynamic) stability and 
suggestions of possible synthesis routes for their formation and subsequent stabilization. The 
molecular species that were examined in detail were chosen from those screened to have a 
predicted theoretical maximum specific impulse in excess of that available from conventional 
chemical systems. It was also intended that this theoretical research effort should be carried out 
with sufficient flexibility to aid in the analysis and guidance of concurrent experimental programs 
of synthesis of high energy mold «lar species. To this end, one task of this research program, 
theoretical studies of the thermochemistry and dissociation reactions of fluorine azide, was 
carried out in parallel with an experimental research study of azide chemistry that was conducted 
at the Rockwell Science Center. In addition, a study of the potential of fluorine azide as a source 
of NFfa 1 A ] for a chemically pumped laser, was carried out for AFWL in conjunction with a 
parallel experimental program at Rockwell Science Center. Finally, a joint theoretical study of 
U3H as a potential fuel was carried out in collaboration with the University of Alabama. 

Theoretical electronic structure calculations were carried out for several high energy species 
that appear to be attractive candidates for advanced chemical propulsion systems. Using 
deliverable specific impulse and storability as the major criteria for the evaluation of new 
oxidizers or fuels, primary consideration was given to ground state molecular structures, of low 
molecular weight, which exhibit a high positive heat of formation. Calculations to date have been 
carried out on: 1 ) light element C3V and C2v structures (H4, U3H, UH3, Li 4 ); 2 ) azide-like 
structures (FN3, a-N202, FNCO, CO3, FNBF); 3 ) cyclic boron structures (B3H3, B2H2NH); 
and 4 ) hypervalent structures (NF5, PF5). 

We find that the ground states of H 4 (tetrahydrogen), UH3 and Li 4 , as C3 V structures, are 
unstable toward dissociation into diatomic fragments. These molecules do not look promising as 
storable, high energy chemical species. We have studied the U3H system in both C3V (trigonal 
pyramid) and C2v (planar kite-like) symmetries. Although both structures are vibrationally stable, 
with the C2v symmetry slightly lower in energy at the MP 2 / 6 - 311 G** level of theory, an analysis of 
the bond character in the C3V and C2 V structures indicate that only the latter is thermodynamical¬ 
ly stable. Using either O2 or F2 as an oxidizer, U3H appears to give about 10 % improvement in 
I sp over LiH. The optimum I sp for the U3H/H2/O2 system is 456 seconds, which can be compared 
with 441 seconds for the UH/H2/O2 system and 441 seconds for the H2/O2 system run at P c = 

800 psia, A/A* = 20 : 1 . 
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All of the azide-like structures are vibrationally stable with high ( 400-500 kJ/mol) positive 
heats of formation. The best candidate appears to be asymmetric dinitrogen dioxide, a-N2C>2, 
which exhibits a heat of formation of + 500 kJ/mol and can deiiver 450 seconds I sp as a hybrid 
propellant with hydrogen, and advanced performance in solid propellant formulations as a 
replacement for AP The energetics and pathways for decomposition of FN3 have also been 
studied. Our calculations, carried out to compliment the experimental studies of FN3 at Rockwell 
Science Center, indicate that NF[a l A ] should be efficiently produced by the thermal decomposi¬ 
tion of FN3. This shows promise as a method to produce a short wavelength chemical laser. The 
structure of this molecule, its predicted vibrational spectrum and the role of low-lying triplet 
states have been examined. 

None of the cyclic boron structures examined to date, which include B3H3 and B2H2NH, are 
highly energetic. Their chemistry is further complicated by condensation reactions which cause 
conversion to higher molecular weight structures. 

The hypervalent NF5 molecule has been studied as a D3h structure, as the ionic NF4 + F" 
form and as an inverted C3V structure. Detailed calculations indicate stability for NF5 only as a 
D 3h structure. Correlated energy calculations were carried out for NF5 and for several azide 
structures. Our analysis indicates that NF5 is a thermodynamically stable compound with a 
predicted heat of formation AHf° (OK) = + 25.1 kJ/mol. An H2/F2/NF5 propellant system is 
predicted to have an optimum I sp of 466 seconds (P c = 800 psia, A/A* = 20 : 1 ), yielding no 
significant improvement over the H2/F2 system. However, if the physical properties of NF5 
indicate the possibility of non-cryogenic storage, NF5 may be a promising candidate for an earth 
storable oxidizer. 









PREFACE 


This report was prepared by the United Technologies Research Center, East Hartford, 
Connecticut, under Contract F04611-8t>-C-00?l. The research was funded by both the Air Force 
Astronautics Laboratory and the Air Force Weapons Laboratory. 

Inclusive dates of research were 1986 August 15 through 1989 October 31. Captain William A. 
Sowell, Lt. Walter J. Lauderdale, Lt. Steven D. Thompson, and Lt. Pete Dolan have served in 
succession as Project Managers for this contract. 

Very useful discussions with Drs. David Benard (Rockwell Science Center), Michael Bowers 
(University of California, Santa Barbara), K. Lammertsma (University of Alabama), D. W. Setser 
(University of Kansas) and Capt. G. P. Perram (AFWL), are also acknowledged. 


All aspects of the research work reported herein were aided by the skilled help of Judith B. 
Addison (UTRC), who carried out much of the computer program development and assisted in the 
analysis of the calculated data and in the preparation of this final report. 



Acoesslon For 


HTIS GRAM 
MIC TAB 
Unannounced 
Justification. 


t 

□ 


Distribution/ 


Availa bility Codes 
Avail and/or 
Special 


Diet 


(d 


V 





TABLE OF CONTENTS 


Section .Page 

SUMMARY. iii 

PREFACE. v 

LIST OF FIGURES . ix 

LIST OF TABLES . x 

1 INTRODUCTION.. I 

2 MATHEMATICAL BACKGROUND . 6 

2.1 Method of Ab Initio Calculation . 7 

2.1.1 Born-Oppenheimer Separation. 7 

2.1.2 Variational Methods. 6 

2.1.3 Configuration Selection. 11 

2.1.4 Multiconfiguration - Self Consistent Field Method (MC-SCF)- 13 

2.1.5 Many-Body Perturbation Theory (MBPT) . 14 

2.2 Ab Initio Gaussian Wavefunction Electronic Structure Codes . 16 

2.2.1 GAUSSIAN 88 . 16 

2.2.2 GAMESS. 18 

2.2.3 CADPAC. 19 

2.2.4 HONDO. 21 

2.2.5 COLUMBUS. 22 

2.3 Spin-Projected Unrestricted Hartree-Fock Method. 23 

2.4 Transition Probabilities . 27 

3 DISCUSSION AND RECOMMENDATIONS. 36 

3.1 Light Element High Symmetry Structures. 44 

3.1.1 H 4 . 44 

3.1.2 Li 3 H . 45 

3.2 Halogen Azides and Azide-Like Structures . 51 

3.2.1 a-N 2 0 2 . 56 

3.2.2 FN 3 . 63 

3.2.3 C1N 3 , BrN 3 . 69 

3.3 Cyclic Boron Structures. 72 

3.4 Hypervalent Structures . 75 

vii 





































TABLE OF CONTENTS (CONTINUED) 


Section 

4 

5 


Bags 


CONCLUDING REMARKS . 82 

REFERENCES. 84 

APPENDIX A - PUBLICATIONS AND PRESENTATIONS. A-l 

APPENDIX B - ABSTRACTS. B-l 

APPENDIX C - THE STRUCTURE AND BONDING OF 

Li 3 H ION-PAIR STATES .... C-l 


viii 












LIST OF FIGURES 


Page 


Eigum 

1 Specific Impulse Obtainable from Fuel/Oxidant 

Energy Release with Hydrogen as a Working Fluid. 39 

2 Low-Lying Quintet States of He 2 . 42 

3 Generalized Reaction Surface . 43 

4 High-Energy Light Element Cj v Molecular Structures. 47 

5 Li 3 H/F 2 Vacuum Specific Impulse . 52 

6 LijH /02 Vacuum Specific Impulse. 53 

7 Energetics of the Chemical Synthesis of Dinitrogen 

Dioxide (a-N20 2 ) . 58 

8 a-N 2 0 2 Dissociation Pathway . 60 

9 H 2 /a-N 20 2 Vacuum Specific Impulse . (2 

10 FNj -* NE + N 2 Reaction Path. 71 

11 Cyclic Boron Compounds. 74 

12 Comparison of D 3 t, and C 2 v Structures of NF 5 . 77 


ix 

















LIST OF TABLES 


I 


Table Pa ge 

1 Ab Initio Structure and Calculated Energy for H 4 

in C 3 V Symmetry . 46 

2 Theoietical Resuits fo - L^H Cj v Structure. 48 

3 Theoretical Results for L^H C& Structure. 50 

4 Ab Initio Structure and Calculated Energy for Li 4 

and LiHj in Cj. Symmetry . 54 

5 Optimized Geometries of Azides. 55 

6 N 2 O 2 Theoretical Predictions . 59 

7 FN 3 Equilibrium Geometry . 65 

8 FN 3 Equilibrium Vibrational Frequencies. 65 

9 FN 3 -* N 2 + NF Transition State . 67 

10 FN 3 -* N 2 + NF TS Vibrational Frequencies. 67 

11 FN 3 -* N 2 + NF Barrier Height . 68 

12 FN 3 3 A” UHF/6-31G* Energy . 68 

13 FN 3 SCF IRC (Reaction Path). 70 

14 FN 3 1 A' - 3 A" Higher Order Energies at SCF Geometry. 70 

15 Ab Initio Structure and Calculated Energy for NF 5 

in D 3h Symmetry 78 

16 Calculated Harmonic Frequencies for NF 5 as a D 31 , Structure. 79 


I 


x 





















SECTION I 


INTRODUCTION 

Conventional chemical systems for rocket propulsion have a predicted theoretical maximum 
specific impulse of approximately 389 seconds (H 2 /O 2 system, 1000 psia - sea level). Although 
some improvements in specific impulse may be possible through development of new chemical 
oxidizers, no chemical system currently under development appears likely to provide a large 
increase above this value. On the other hand, advanced space missions presently considered need 
the development of propulsion systems with specific impulses of 800 to 2,000 seconds for earth 
orbiter travel, and much higher specific impulses for interplanetary travel. Actual operation of 
such systems clearly remains in the future, but exploratory studies of such systems certainly are 
needed at the present time. The search for high-performance rocket propellants has, therefore, of 
necessity turned to consideration of advanced chemical systems, including various metastable 
species such as He/ and He 2 *. compounds with large positive heats of formation (such as FN 3 ), 
and antimatter propulsion. Included among these species are stabilized free radicals, which have 
received considerable attention in the past. Performance considerations always make light element 
molecules the most promising species for advanced rocket propulsion systems and as a result, a 
number of efforts have been undertaken to characterize molecular structures of light element 
species that have a high energy content because of the peculiar nature of their chemical bonding. 

It was the goal of this research program to identify those light element molecular structures with 
high potential as advanced propellants and to characterize their energy content, possibl *'nthesis 
routes and their chemical stability. 

The advantages of using highly energetic compounds or radicals as propellants for rocket 
engines have been discussed by several authors (References 1-3), and the performance capabilities 
of such materials as the light metal hydrides or atomic hydrogen are reasonably well established. 
In principle, the stabilization of atomic species appears simple, and depends on keeping the 
atoms physically separated, thus preventing them from recombining into their thermodynamically 
stable molecular state. In practice, this may be achieved by dispersing the atoms in a solid matrix 
at temperatures low enough to prevent their diffusion through the matrix. Considering isolation of 
the atoms as the criteria, one could expect stabilization of up to 50% atoms for a perfect body- 
centered cubic arrangement. The manner of preparation by condensation from a gaseous 
discharge, however leads to a disordered random configuration, in which one might expect 
concentrations of approximately 10% (References 4 and 5). Other calculations have indicated an 
upper limit of about 3% (Reference 6 ). Atomic hydrogen as well as other atoms have been 
stabilized in a solid matrix at liquid helium temperatures (4 K) (References 6-9). The observed 
concentrations have been smaller than predicted by theory and are generally of the order of 0.2 to 
0 . 001 %. 
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At temperatures above 4 K. most studies have shown that diffusion sets in rapidly. Recombi¬ 
nation then occurs and is generally complete at 20-30 K. This behavior is as expected considering 
the small size of the atoms and the nature of the matrices. Hydrogen atoms would therefore not 
be expected to have any long-term stability at temperatures near 77 K (the normal boiling point of 
liquid nitrogen). Similar behavior has been observed for other light element free radicals (Refer¬ 
ence 10). In summary, efforts at stabilization by physical methods, such as matrix isolation, of 
significant quantities of dissociated atoms, in order to take advantage of their large energy of 
recombination have been quite discouraging, as evidenced by both theory and experiment. 

The potentialities of another class of possible compounds in which helium is one of the 
constituent atoms has also been well explored. There are two types of compounds that might be 
formed from the noble gases such as helium. The first type is the so-called van der Waals 
molecule which forms when an inert atom such as helium combines with another atom or 
molecule. In this case the separate atoms, which are ail in their ground states, form complex 
structures due to weak van der Waals attractive forces. It has been shown (Reference 11) that the 
hydrates, solvates, and mercury compounds of helium are loosely bound molecular complexes held 
together by forces of the weak van der Waals type. Complexes of this type have small heats of 
dissociation ( < 2.0 kjoule/mol) and are of little use as potential energy sources. 

The second type of compound that might be formed from the noble gases is the metastable 
structure that results when an inert atom is electronically excited by an electric arc (or other 
means) and combines with another atom or molecule. The possibility of forming compounds in 
this way arises from the supposition that, although the normal helium atom is chemically inert, 
helium excited to higher electronic states might behave as an alkali metal. The energy available in 
electronically excited states is often quite large (for example, the lowest excited state of helium, a 
triplet state, is approximately 1900 kjoule/mol above the ground state) and an atom excited to one 
of these levels and stabilized by compound formation might retain a substantial amount of this 
energy for subsequent release in a controlled reaction such as a chemical rocket. A most attrac¬ 
tive molecule for such an application would be the metastable a 3 Z u + state of the helium 
molecule since the molecular weight of the decomposition products is four and the average value 
of y approaches 1.67, thus insuring efficient energy conversion in the exhaust nozzle. 

Stabilization of compounds of excited helium is not likely at ordinary temperatures, but it 
has been suggested that low temperature stabilization might prove feasible if the average lifetime 
of the excited molecule is sufficient to allow trapping of the energy in a solid matrix or by other 
means. In order to determine the feasibility of storing an energy-rich molecule, an exact analysis 
of the energy levels and expected lifetimes of the excited states is required. This Research Center 
has carried out such a study of metastable helium compounds for the Air Force Rocket Propul¬ 
sion Laboratory (Reference 12) in the early sixties. Detailed quantum mechanical studies of 
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several metastable helium compounds (Hei. HeH, HeBe) were carried out to determine the energy 
levels and lifetimes of these species. These early calculations were at the then state-of-the-art and 
included some of the very first configuration - interaction (Cl) studies ever carried out on excited 
states of diatomic molecules. After a lengthy analysis of the deexcitation mechanisms of such 
metastable species, it was concluded that stabilization was not likely, even in matrices at 4 K, and 
efforts to synthesize such compounds were phased out. 

The idea of stabilizing metastable helium atoms by formation of a spin-polarized solid 
phase has also been proposed (References 13 and 14). This idea apparently developed from a 
suggestion by Garrison, et al. (Reference 15) that the lowest 5 2 g + state of He 2 , formed from two 
metastable He atoms, exhibited weak chemical binding (0.07 eV) and no first order radiative 
decay process to the underlying continuum states. Arguments were given (Reference 13) that such 
a material should form an insulating classical solid and that the metastability could be trapped by 
either intrinsic effects of the condensed phase or by external influences such as pressure or 
electric or magnetic fields. 

In an effort to analyze these proposed spin-polarized He compounds, this Center carried 
out a detailed quantum mechanical study of molecule formation through collisions of triplet 
metastable helium atoms (Reference 16). Contrary to a previous theoretical study (Reference 15), 
the lowest lying molecular state was found to be of 5 fl g symmetry and all of the low-lying quintet 
states of He 2 were shown to arise from the ionic He" - He + interaction. The 0 g + and lg 
components mix via spin-orbit coupling to the autoionizing singlet and triplet manifolds of He 2 
with a cross-sectior of ~ 100 A 2 , The 0 g ", 2g and 3g components mix via spin-orbit coupling to 
the ionic He" - He* manifold with radiation lifetimes of ~ 10 tx sec to the He 2 + + e underly¬ 
ing continuum. Thus, again we found no evidence that metastable helium compounds could be 
stabilized against radiative decay or autoionization processes. 

The past emphasis on free radical stabilization in atomic species or metastability in helium 
compounds has caused many investigators to overlook the possibility of finding high-energy 
species, possibly suitable for propellant application, which represent molecules with steric atomic 
arrangements that represent localized minima on the ground state potential surface. To a certain 
extent, this has been due to the enormous effort required in any experimental program of 
screening the large number of possible molecular conformations which can be made up from, for 
example, the lightest four or five elements. However, theoretical techniques do exist for screening 
this large aggregate of possible molecules to determine those which might be sufficiently stable 
and energetic for propellant applications. More definitive calculations could then be made for the 
most promising of these species to verify the results of the screening process and to provide 
additional data on which to base recommendations for experimental synthesis studies. 
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In view of the above considerations, and because of the high performance potential of such 
systems, a broad, comprehensive research investigation of light-element high-energy compounds 
that might be useful for advanced chemical propulsion systems was undertaken by this Center for 
AFRPL in the late sixties (Reference 17). In particular, this study consisted of a theoretical 
evaluation of the possibility of high-energy propellant formation from each of the diatomic and 
symmetric triatomic species which could be formed from elements having atomic numbers of 
eighteen or less. The objective of the study was the determination of those stable or metastable 
chemical systems, included in the above categories, which are most likely to provide specific 
impulses of the order of 500 seconds or greater, and, as such, are most worthy of further theoreti¬ 
cal and experimental investigation. Initial emphasis was mainly upon highly-excited metastable 
states, theoretically capable of delivering very high impulses. During the latter stages of this work, 
the study was broadened to concentrate more heavily on energetic ground states which might 
possibly be more easy to stabilize. Over 10.000 electronic states of nearly 500 molecular species 
were screened in that study 

In the first phase of that study emphasis was placed upon the possible use of the light, 
noble-gas difluorides as advanced oxidizers. In this phase, calculations were performed for the 
low-lying states of HeF 2 and NeF 2 and their diatomic and atomic fragments. All states of the 
difluorides which were investigated were found to be repulsive, as were the states of HeF and 
NeF. It was concluded that the light noble-gas fluorides do not constitute an attractive source of 
high-energy oxidizers. In the second phase of the investigation, emphasis was placed upon 
ascertaining the stability of a number of compounds believed to be metastable or to have a high 
positive heat of formation. Calculations were performed for the ground states of the polyatomics 
FI^Ne, H 2 Ne, F^Na, F^Ar, C^Ne, various possible diatomic and atomic fragments thereof, and 
for several low-lying states of the diatomics ONe and OAr. Except for the hydrides, such as 
H 2 Ne, no appreciable bonding was exhibited for any molecular arrangements of the ground state 
of these species. It was concluded, at the time, that although several of these hydrides showed 
promise as high-energy species, most were ruled out for their potential use as fuels owing to their 
molecular weight. 

Very recently, there has been a renewed interest in light element compounds in geometric 
arrangements other than the stable ground state form. Several authors (References 18-22) have 
discussed various forms for the F4 molecule, sometimes referred to as tetrahydrogen. The systems 
HeFl 2 and Li m H n have also been studied both experimentally and theoretically (Refer¬ 
ences 23-30). In all of these systems, molecular arrangements have been found which correspond 
to an energetic local minima on the ground state potential energy surface of the reactant atoms. 
These high energy molecular arrangements can be referred to as strained isomeric forms. In many 
cases, these high energy isomeric species exhibit ionic binding, whereas the lowest energy 
arrangement is a covalent structure. There are many such isomeric arrangements known in 
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organic chemistry, the vinylidine form of acetylene, being a typical example. A survey of strained 
organic molecules has been given by Newton (Reference 31). Under the present research program 
we have carried out a systematic study of such isomeric arrangements of light element com¬ 
pounds, using the most advanced quantum mechanical techniques that are available, to screen 
those having the highest potential as advanced propellant types and to suggest possible reaction 
routes for their synthesis and stabilization. The results of our preliminary studies of such 
compounds and their performance evaluation are given below in the Discussion and Recommen¬ 
dations section. 

The general composition of this report is as follows. In Section Z we present a description of 
the mathematical methods which were employed in this research. Included in Section 2 are 
sub-sections which deal with the construction of electronic wavefunctions, the calculations of 
expectation properties, the evaluation of molecular transition probabilities, and the calculation of 
electronic wavefunctions using ab initio methods. The calculated results and pertinent discussion 
are presented in Section 3. Based on these present studies, we have identified several areas for 
further technical work. Recommendations are also presented in Section 3 and a concluding 
summary of this program is given in Section 4. A list of all technical presentations and publica¬ 
tions that have resulted from this research program are given in Appendix A. Abstracts of the 
pertinent published papers are then collected in Appendix B. Appendix C contains a recent 
publication on L^H. 
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SECTION 2 


MATHEMATICAL BACKGROUND 

Central to these theoretical studies are the actual quantum-mechanical calculations which 
must be carried out for the atomic and molecular species. In particular, the methodology for 
computing the minimum energy reaction pathways and testing for vibrational stability must be 
carefully analyzed. For added clarity, various aspects of these calculations are discussed below. 

Much evidence on diatomic and polyatomic systems indicates the inadequacy of minimum 
basis sets for constructing quantitatively correct molecular wavefunctions (References 32 and 33). 
This means inner-shell and valence-shell orbitals of quantum numbers appropriate to the atoms 
(Is. 2s, 2p, for C, N, O; etc.). The main deficiency of the minimum basis set is its inability to 
properly describe polarization and the change of orbital shape for systems which exhibit large 
charge transfer effects. Values of the screening parameters for each orbital can either be set from 
atomic studies or optimized in the molecule; the latter approach is indicated for studies of higher 
precision. When high chemical accuracy is required, as for detailed studies of the ground or a 
particular excited state of a system, a more extended basis must be used. Double-zeta plus 
polarization functions or optimized MO’s are usually required as a minimum representation for 
reliable calculated results of chemical accuracy. 

The chosen basis sets give good results only when used in a maximally flexible manner. This 
implies the construction of perturbation expansions or the use of Cl wavefunctions with all kinds 
of possible orbital occupancies, so that the correlation of electrons into overall states can adjust to 
an optimum form at each geometrical conformation and for each state. Except when well-defined 
pairings exist, as for closed shell and exchange dominated systems, a single-configuration study 
(even of Hartree-Fock quality) will be inadequate. 

Proper electronic states for systems composed of light atoms should possess definite 
eigenvalues of the spin operator S 2 as well as an appropriate geometrical symmetry. The geometri¬ 
cal symmetry can be controlled by the assignment of orbitals to each configuration, but the spin 
state must be obtained by a constructive or projective technique. Formulas have been developed 
(Reference 34) for projected construction of spin states from orthogonal orbitals, and programs 
implementing these formulas have been in routine use at UTRC for several years. One of the least 
widely appreciated aspects of the spin-projection problem is that the same set of occupied spatial 
orbitals can sometimes be coupled to give more than one overall state of a given spin quantum 
number. It is necessary to include in calculations all such spin couplings, as the optimum coupling 
will continuously change with changes in the molecular conformation. This is especially important 
in describing degenerate or near-degenerate excited electronic states. 
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In the sections below, we describe the several mathematical approaches that are applicable 
to calculation of the electronic structure of molecules, to the calculation of potential energy 
surfaces for chemically reacting systems, and to the subsequent calculation of radiation and 
collisional processes. Since several different approaches are indicated, we describe their expected 
areas of applicability. 

2.1 Method of Ab Initio Calculation 

2.1.1 Born-Oppenheimer Separation 

For a system of n electrons and N nuclei, and considering only electrostatic interactions 
between the particles, we have for the total Hamiltonian 


V r>2 ^ 2 

*-*■<-25^+5; 


X X v.-*> + 2 jr i vv, + XI v ' v i 


L 1 a= l 


a = 1 i = 1 


.= lj=l 
« 


( 1 ) 


where 


K.I = -^X V ? +V "(f..RN) 


( 2 ) 


and where m e , m a , Mj, are the masses of the electron, atom and combined system mass, respec¬ 
tively. Now since the ratios m e /m a , and m e /M T are both small, (2 • 10~ 6 - 5 • 10~ 4 ) we can effect 
a separation of the electronic and nuclear coordinates treating the total wavefunction as a product 
of a nuclear and an electronic part. We have 


V<r„,R N ) = XZk(RN>k(fn,RN) 

k 


(3) 


where ^{x n , Rn) is an electronic wavefunction parametric in the nuclear coordinates as given in 
Equation (3) and &(R N ) are nuclear motion wavefunctions which satisfy (neglecting terms of the 
order of m e l m a ) 
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(4) 



The cross term in V a • can be eliminated by a proper change of variables and Equation 
(4) then reduces to a 3N-3 dimensional Schrodinger equation. 

For most systems, where the velocity of motion of the nuclei is slow relative to the electron 
velocity, this decoupling of electron and nuclear motion is valid and is referred to as the adiabatic 

approximation. Equation (3) thus defines an electronic eigenstate »/\(r„, Rn), parametric in the 

nuclear coordinates, and a corresponding eigenvalue Ei^Rn) which is taken to represent the 
potential energy curve or surface corresponding to state k. 

2.1.2 Variational Methods 

By an ab initio method is meant one that starts from a zero-order Hamiltonian which is 
exact except for relativistic and magnetic effects, and which involves the evaluation of electronic 
energies and other relevant quantities for wavefunctions which are properly antisymmetrized in 
the coordinates of all the electrons. For a system containing n electrons and M nuclei, the 
zero-order Hamiltonian depends parametrically on the nuclear positions and is of the form 


= 


1 N n M M . 

I AV + X A 

i = l i = lj=l Tj - Rj lsicj Rj - Rj 1 SI<J | r i r j| 


(5) 


where Z, and R, are the charge and position of nucleus i, r, is the position of electron j, and Vj 2 

is the Laplacian operator for electron j. All quantities are in atomic units, i.e. lengths in bohrs. 
energies in hartrees (1 hartree = 2 Rydbergs). The many-electron wavefunction consists of one, or 
a linear combination, ip = Z c^, of terms of the form 


n 

%(R> = ja n*RR&M 

i = 1 


( 6 ) 


where each (p^ is a spatial orbital, A is the antisymmetrizing operator, O s is the spin-projection 
operator for spin quantum number S, and 9,, m is a product of a and ft one-electron spin 


8 




functions of magnetic quantum number M s . No requirement is imposed as to the double occu¬ 
pancy of the spatial orbital, so linear combinations of the form given by Equation (6) can describe 
a completely general wavefunction. The spatial orbitals <p^ may be whatever basis orbitals have 

been introduced, arbitrary linear combinations thereof, or specific linear combinations deter¬ 
mined pursuant to the particular calculational method in use. 

The spatial orbitals <f> M i , the spin functions 6 and the coefficients of different xp ^, if a 
linear combination of is used, may be explicitly determined by invoking the variational 

principle. Various specific methods are described below for determining wavefunctions. However, 
we should first observe that the adequacy of an ab initio calculation, or for that matter any energy 
calculation, will depend crucially upon the extent to which the wavefunction can be qualitatively 
appropriate. Some of the considerations surrounding the choice of a wavefunction are the 
following: 

(i) necessity that the wavefunction possess sufficient flexibility to be able to describe 
dissociation to the correct atomic and molecular fragments as various internuclear separations are 
increased; 

(ii) maintenance of equivalent quality of calculation for nuclear geometries differing in the 
nature or number of chemical bonds; 

(iii) ability to describe degenerate or near-degenerate electronic states when they are 
pertinent; 

(iv) ability to describe different electronic states to equivalent accuracy when their interrela¬ 
tion (e.g., crossing) is relevant, in particular, ability to describe ionic-valence state mixing; 

(v) ability to represent changes in the coupling of electron spins as bonds are broken or 
reformed. 

The foregoing considerations indicate that it will often be necessary to consider wavefunc¬ 
tions with more than a minimum number of singly-occupied spatial orbitals, and that there will 
be many potential curves or surfaces for which a wavefunction consisting of a single ip M cannot 

suffice. It will then be necessary to allow mixing of % with different degrees of orbital spatial 

occupancy so as to obtain smooth transitions from the occupancies characteristic of separated 
atoms or molecules (or ions) to those characteristic of a compound system or a different fragmen¬ 
tation. 

Another implication of the considerations surrounding the choice of a wavefunction is 
related to the treatment of electron spin. Not only is it necessary to require that the wavefunction 


9 




be an eigenfunction of S 2 and S z but it is also necessary to take account of the fact that under 
many conditions, there will be more than one spin eigenfunction of given S and m s . The different 
spin eigenfunctions correspond to different couplings among the individual spins. Since reactive 
processes involve the breaking and forming of electron-pair bonds, they must necessarily be 
accompanied by reorganizations of the spin coupling. A failure to take account of this will lead to 
qualitatively inappropriate wavefunctions. 

In Hartree-Fock calculations ip (R) is restricted to a single ip M which is assumed to consist 
as nearly as possible of doubly-occupied orbitals. The orbitals (p^ are then selected to be the 
linear combinations of basis orbitals best satisfying the variational principle. Writing 
<pn j = 2 a y# v , the a M are determined by solving the matrix Hartree-Fock equation 


X F ^v, = €j X < eaCh *) ( 7 ) 

V V 


where t, is the orbital energy of <p^ . 

The Fock operator F iv has been thoroughly discussed in the literature (Reference 35) and de¬ 
pends upon one- and two-electron molecular integrals and upon the a* . This makes Equa¬ 
tion (7) nonlinear and it is therefore solved iteratively. UTRC has developed programs for solving 
Equation (7) for both closed and open-shell systems, using basis sets consisting of either Slater- 
type or Gaussian-type atomic orbitals. Examples of their use are in the literature (Reference 36). 

In configuration interaction calculations, the overall wavefunction has more than one term, 
tpf,, and the c^ are determined by invoking the variational principle to obtain the secular 

equation 


X(H„v-ws„ v )c v = 0 (each pi) 

V 

where 


( 8 ) 


H u 


-h 


(R^R^RVr 




(R)^R>/r 


(9) 
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Equation (8) is solved by matrix diagonaiization using either a modified Givens method (Refer¬ 
ence 37) or a method due to Shavitt (Reference 38) or Raffenetti (Reference 39). 

The matrix elements H^„ and may be reduced by appropriate operator algebra to the 
forms 


■i n 

H„v = 5> p (0m|O,P|0m)( F1%>( f -Rn)|k(Rn)P| n ,y "(G^N)} GO) 


i = 1 


V = 2 f p (tfM|0 1 P|d M )(n%i(Fi,R N )|P| fl^(F„R N )) ( 11 ) 

p .=i i=i 

where P is a permutation and e p its parity. The sum is over all permutations. (0^O % P\0 m ) is a 
“Sanibel coefficient” and the remaining factors are spatial integrals which can be factored into 
one- and two-electron integrals. If the are orthonormal. Equations (10) and (11) become more 
tractable and the H^ v and S^„ may be evaluated by explicit methods given in the literature (Ref¬ 
erence 34). Computer programs have been developed for carrying out this procedure, and they 
have been used for problems containing up to 106 total electrons, 10 unpaired electrons, and sev¬ 
eral thousand configurations. 

The Cl studies described above can be carried out for any orthonormal set of <p m for which 
the molecular integrals can be calculated. Programs developed by UTRC make specific provision 
for the choice of the <pp, as Slater-type atomic orbitals, as Gaussian-type orbitals, as symmetry 

molecular orbitals, as Hartree-Fock orbitals, or as more arbitrary combinations of atomic 
orbitals. 

The one- and two-electron integrals needed for the above described method of calculation 
are evaluated for STO’s by methods developed by this Center (Reference 40). For Gaussian 
orbitals, either the Camegie-Mellon integral package (Reference 41) or the integral routines 
incorporated in the HONDO or GAMESS programs (Reference 42) can be employed. 

2.1.3 Configuration Selection 

Using a double-zeta plus polarization basis set of one-electron functions, a typical system 
can easily have of the order of 10 6 configurations in full Cl (that resulting from all possible orbital 
occupancies). It is therefore essential to identify and use the configurations describing the 
significant part of the wavefunction. There are several ways to accomplish this objective. First, 
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one may screen atomic-orbital occupancies to eliminate configurations with excessive formal 
charge. Alternatively, in a molecular-orbital framework, one may eliminate configurations with 
excessive numbers of anti-bonding orbitals. A third possibility is to carry out an initial screening 
of configurations, rejecting those whose diagonal energies and interaction matrix elements do not 
satisfy energy significance criteria. 

Another common method of classifying configurations is to examine the total number of 
orbitals in the wavefunction that differ from the SCF reference wavefunction. We write the Cl 
wavefunction as 


ipc i = C o 0o + + X Cj$j D + ••• 

' j 


( 12 ) 


where a single excitation function, differs from the SCF reference by one orbital and the 
double excitation function^ by two, and so on. The Cl coefficients, C„ are then determined 

variationally to yield the lowest possible total energy. In the limit of a complete basis set (N—oo), 
and where all possible substitutions are included in . the variational energy approaches the 
correct nonrelativistic Born-Oppenheimer result Errors arise from a truncation of the functions 
used to determine the SCF reference wavefunction and from the truncation of the excitation func¬ 
tions series in V'ci - The reference wavefunction \p Q will typically be the same for a Cl or Many 
Body Perturbation Theory (MBPT) calculation; however, for a molecule even as small as water. 
ipci becomes a function of a very large number of basis functions. Additionally, the reference 
state itself may be multi-dimensional in order to describe systems such as diradicals or systems 
with orbital degeneracy. Because of this, one must truncate the expansion and eliminate unimpor¬ 
tant configurations. In general one can show that 


(VoW <p J ) = (vM>°) = o 


(13) 


In words, the matrix elements between the reference wavefunction and triple and quadruple exci¬ 
tations is zero. Thus to first order, only single and double excitations contribute. Because of this, 
many Cl calculations attempt to include all single and double excitations in the expressions for 
ipa ■ To go beyond this, generally more than one reference wavefunction is used. Programs to han¬ 
dle configurations on all the above criteria are available at UTRC. 
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Other, potentially more elegant methods of configuration choice involve formal app.oaches 
based on natural orbital (Reference 43) or multiconfiguration SCF (Reference 44) concepts. To 
implement the natural-orbital approach, an initial limited Cl wavefunction is transformed to 
natural orbital form, and the resulting natural orbitals are used to form a new Cl. The desired 
result is a concentration of the bulk of the Cl wavefunction into a smaller number of significant 
terms. The multiconfiguration SCF approach is more cumbersome, but in principle more 
effective. It yields the optimum orbital choice for a preselected set of configurations. This 
approach works well when a small number of dominant configurations can be readily identified. 
The method is described briefly below. 

2.1.4 Multiconfiguration - Self Consistent Field Method (MC-SCF) 

The Hartree-Fock self consistent field method has been proven to be a powerful tool for the 
calculation and understanding of many ground state properties of molecules in the vicinity of 
their equilibrium structure. However, in most cases the one configuration Hartree-Fock approxi¬ 
mation is not adequate to properly describe the dissociation of molecular bonds. Also, many 
excited states cannot be represented by a single configuration wave function. In order to calculate 
properties for such states, or to investigate the formation of molecular bonds, one often needs 
multiconfiguration wave functions for which both the linear coefficients, C, of the configuration 
expansion 


V ~ X C 'W 


(14) 


as well as the set of orthonormal molecular orbitals {^j} , from which the configurations are 
constructed, are optimized according to the variational principle. As is well known, this "MC- 
SCF” problem presents many more difficulties than the simple one configuration Hartree-Fock 
case and much work has been devoted to obtaining convergent solutions during the last decade. 

The difficulties mainly arise from the fact that for general MC-SCF wave functions the 
energy is not invariant with respect to rotations between occupied orbitals. Hence, instead of a 
relatively simple pseudo-eigenvalue equation in the one determinant case, the set c r coupled Fock 
equations 


X F 'M) = (15) 

i i 
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with the hermiticity conditions 


has to be solved. e,j are Lagrange multipliers which account for the orthonormality constraints 
imposed on the orbitals. The Fock operators F,j depend on the orbitals (0[) and the set of Cl co¬ 
efficients (CjJ. In analogy to the one determinant case many attempts have bee., made to solve 

these equations iteratively by keeping the Fock operators fixed in each iteration step. Then the 
Lagrange multipliers can be expressed by coupling operators constructed such that the Fock 
equations are transformed into pseudo-eigenvalue equations yielding the improved orbitals. These 
are used in a second step to determine new Cl coefficients by diagonalizing the Cl matrix 

. The convergence of these algorithms, however, has often been found to be poor. 


A second group of MC-SCF methods is based on the generalized Brillouin theorem. In these 
methods, the orbital changes are derived from the coefficients of a Cl expansion consisting of the 
MC-SCF wavefunction and ail one-electron singly excited configurations. A computer program 
(ALIS) implementing this method has been developed by Ruedenberg, et al. (Reference 45). A 
somewhat more elegant program (GAMESS) which also incorporates analytical gradients is also 
available (Reference 42). More recently, various methods have been proposed which are based on 
direct minimization of the energy, avoiding the Fock operators altogether. Such methods are now 
being studied in our laboratory and will be incorporated into our existing computer programs if 
they prove to be highly efficient. 

2.1.5 Many-Body Perturbation Theory (MBPT) 

In MBPT we again begin with the SCF wavefunction (ip Q ) as our reference and attempt to 
account for E c , the correlation energy. The concept of excitation functions described in the above 
section on Cl calculations carries over to MBPT calculations. In MBPT one can write the 
wavefunction, ip p , as 


</>P = e 7 ^} 


(17) 
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where T is an excitation operator defined as 


T = Ts + T d + Tt + Tq... (18) 

where S, D, T. and Q refer to single, double, triple, and quadruple substitutions respectively. One 
can write T n , in general, where n refers to the number of excitations as 



abc... 


abc 

ijk.. 




(19) 


where a.b.c,... are excited orbitals and i,j,k,... are orbitals occupied in <p 0 . The total energy is now 
given by 


Embpt - (<Po\W>e r \<Po) 


( 20 ) 


To evaluate Embpt. the tf ' from above must be determined. An equivalent expression that 
makes the perturbation expansion clearer is 


Embpt = ^ ^o|K[(E o -Ko) _1 3t] k |0o) n\) 

k =0 

where the sum is over only so-called linked diagrams, Dfco has eigenfunctions <p 0 , and the expan¬ 
sion is of orders in the perturbation, V = (K - 3t„). The k = 0 term gives the reference energy 
and for k > 0 correlation corrections are included. In practice the MBPT total energy is calcu¬ 
lated by truncating the T operator expansion and projecting Ke T "|0 o ) onto the appropriate n- 

space. This leads to a set of nonlinear coupled equations for the t* coefficients which 

correspond to the Cl expansion coefficients. The equations are solved iteratively and E M bpt 
evaluated. In practice T 4 is an upper limit that corresponds to quadruple substitutions. The series 
is an oscillatory convergent sum, which in practice has proven to be at least as accurate as E^j 
with single and double excitations included. 

The best possible method to use would be a full Cl (all possible configurations) with a 
complete basis set. However, since the number of configurations is proportional to (nf , where n 
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is the number of basis functions and ( is the level of excitation, it would be prohibitive to even 
include all single and double excitations. This truncation causes the loss of size consistency, which 
implies that the energy calculated for A and B as a molecular system, but dimensionally far apart, 
is the sum of the energy calculated for A and B separately. In a size extensive calculation the 
energy is proportional to the size of the system. These properties are very important if one wishes 
to compute correct relative energies on a potential energy surface, a necessary criteria for defining 
the reactive pathways of interest in this research program. MBPT, on the other hand, is guaran¬ 
teed to have the correct size-dependence because the expansions contain only the so-called linked 
diagrams. In addition, because of the computational efficiency of MBPT calculations can be 
performed up to fourth order in the perturbation expansion and include single, double, triple, and 
quadruple excitations in the calculation of the correlation energy. Such calculations are usually 
performed with at least a split valence plus polarization basis set. For the light element com¬ 
pounds which exhibit ionic bonding, the inclusion of diffuse basis functions and possibly higher 
polarization d or f-functions may be required for a more quantitative treatment. 

2.2 Ab Initio Gaussian Wavefunction Electronic Structure Codes 

Owing to the complexity of evaluating multicenter electron repulsion integrals over Slater- 
type (exponential) orbitals, various groups have adopted a computational approach to electronic 
structure calculations based on gaussian orbitals. A highly developed series of programs, named 
GAUSSIAN 8X (X = 2, 6, 8), is available from Carnegie-Mellon University (Reference 41). A 
second series of programs has evolved from the original version of Dupuis and King's (Refer¬ 
ence 46) HONDO code. This code has been further developed by separate groups as GAMESS 
(Reference 42), CADPAC (References 47) and HONDO 7 (Reference 48). Finally a new code 
named COLUMBUS (Reference 49 and 50) has been developed by the Ohio State/Argonne/Bat- 
telle group. UTRC has been designated as a beta test site for this new development. A brief 
description of the features of these codes follows. 

2.2.1 GAUSSIAN 88 

GAUSSIAN 88 is a connected system of programs for performing ab initio molecular orbital 
(MO) calculations. It represents further development of the GAUSSIAN 70/76/80/82/86 systems 
already published. The contributors to this program include: M. J. Frisch, J. S. Binkley. H. B. 
Schlegel, K. Raghavachari, C. F. Melius, R. L. Martin, J J. P. Stewart, F. W. Bobrowicz. C. M. 
Rohlfing, L. R. Kahn, D. J. Defrees, R. Seeger, R. A. Whiteside, D. J. Fox, E. M. Fleuder. and J. 

A. Pople. GAUSSIAN 8X was originally implemented on the chemistry department DEC VAX 
11/780 computer at Carnegie-Mellon University. Since then this program has been installed on a 
number of different computers. 

GAUSSIAN 88 was designed with a transparent input data stream, making this program 
very user friendly. All of the standard input is free-format and mnemonic. Reasonable defaults 
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for input data have been provided, and the output is intended to be self-explanatory. Mechanisms 
are available for the sophisticated user to override defaults or interface their own code to the 
GAUSSIAN system. In this respect, we intend to utilize GAUSSIAN 88 as a fundamental 
framework for several applications. Options are being incorporated into this code to provide 
capabilities beyond Hartree-Fock and various perturbation theory options. 

The capabilities of the GAUSSIAN 88 system include: 

a) Calculation of one- and two-electron integrals over s, p, d, and f contracted gaussian 
functions. The basis functions can either be cartesian gaussians or pure angular mo¬ 
mentum functions and a variety of basis sets are stored in the program and can be 
requested by name. 

b) Self-consistent field calculations for restricted closed-shell (RHF), unrestricted open- 
shell (UHF), and open-shell restricted (ROHF) Hartree-Fock wavefunctions as well as 
those types of multiconfigurational wavefunctions that fall within the Generalized Va¬ 
lence Bond-Perfect Pairing (GVB-PP) formalism. 

c) Evaluation of various one-electron properties of the Hartree-Fock wavefunction, in¬ 
cluding Mulliken population analysis, multipole moments, and electrostatic fields. 

d) Automated geometry optimization to either minima or saddle points, and analytical or 
numerical differentiation to produce force constants, polarizabilities, and dipole de¬ 
rivatives. This feature can be used to develop minimum energy reaction paths along a 
complicated many dimensional potential energy surface. 

e) Correlation energy calculations using M4>ller-Plesset perturbation theory carried to 
second, third, or fourth order. 

f) Correlation energy calculations using configuration interaction (Cl) with either all dou¬ 
ble excitations (CID) or all single and double excitations (CISD). 

g) Correlation energy calculations using coupled cluster theory with double substitutions 
(CCD). 

h) Correlation energy calculations using quadratic convergence SCF (QCSCF). This is a 
new highly efficient size-consistent method recently developed by Pople. 

i) Analytic computation of the nuclear coordinate gradient of the RHF, UHF, ROHF, 
GVB-PP. MP2, CID and RCISD energies. 

j) Computation of force constants (nuclear coordinate second derivatives), 
polarizabilities, hyperpolarizabilities, dipole derivatives, and polarizability derivatives 
either analytically or numerically. 
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k) Harmonic vibrational analysis. 

l) Determination of intensities for vibrational transitions at the HF, MP2, and Cl levels. 

m) Testing the SCF wavefunctions for stability under release of constraints. 

n) Correlated electron densities and properties. 

o) Minimum-energy pathway following from products to reactants through a transition 
state-intrinsic reaction coordinate finder. 

2.2.2 GAMESS 

A wide range of quantum chemical computations are possible using GAMESS (Refer¬ 
ence 42), a refinement by M. Schmidt and S. Elbert of the original HONDO code (Reference 46). 

The capabilities of this code include: 

a) Calculations of RHF/UHF/ROHF/GVB-SCF molecular wavefunctions. 

b) Calculations of multiconfigi'idtion SCF (MCSCF) wavefunctions. 

c) Calculations of Cl wavefunctions using the unitary group method. 

d) Optimization of molecular geometries using an energy gradient in terms of Cartesian 
or internal coordinates. 

e) Searches for potential energy surface saddle points. 

f) Tracing the intrinsic reaction path from a saddle point to reactants or products. 

g) Computation of normal modes and vibrational frequencies. 

h) Calculation of the following properties: 

1. dipole, quadrupole, and octupole moments 

2. electrostatic potentials 

3. electric field and electric field gradients 

4. electron density and spin density 

5. Mulliken and Lowdin population analysis 

6. localized orbitals by the Boys method 

7. virial theorem and energy components. 
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GAMESS is a synthesis, with many major modifications, of several programs. A large part of the 
program is from HONDO. For pure sp basis sets, the HONDO symmetry and supermatrix proce¬ 
dure has been adapted from GAUSSIAN 76 and GRADSCF integrals, both for the SCF and gra¬ 
dient parts. The GVB section is a heavily modified version of GVBONE. A Boys localization 
algorithm is implemented from a heavily modified version of Streitweiser’s QCPE program. 

The Cl module is based on Brooks and Schaefer’s unitary group program which was 
modified to run within GAMESS, using a Davidson eigenvector method written by S. T. Elbert. 
The MCSCF module is a Newton-Raphson procedure, developed at NRCC, based on the unitary 
group Cl package. The intrinsic reaction coordinate pathfinder was written at North Dakota S?°te 
University. 

2.2.3 CADPAC 

The Cambridge Analytic Derivatives Package (CADPAC) (Reference 47) is a group of 
programs which has been under development at Cambridge University, UK. It originated as a 
version of Dupuis and King’s HONDO program. From its initial state as an SCF gradient 
package for closed-shell and UHF wavefunctions, this program has been extensively modified 
with many of the old features being enhanced and many new features being added. The input data 
is now in free format with a ‘keyword’ system, to make the program easier to use. The integral 
routines use essentially the same methods as those in HONDO, being based on the Rys polyno¬ 
mial method, but have been extended to cover f-functions. These routines have also been 
vectorized to take advantage of the much greater availability of supercomputers such as the Cray 
in recent years. The SCF programs are by now a blend of techniques, but still contain a few 
features from the initial HONDO program, particularly in the way symmetry is handled. There 
have been many modifications to improve efficiency and the addition of level-shifting and 
damping techniques and the implementation of the DIIS method to aid convergence. A high-spin 
open-shell SCF program, and a completely general open-shell SCF program have also been 
included. 

The gradient routines have been considerably altered from those in the original HONDO, 
the method having been changed to one which is more efficient and easier to extend to higher-or¬ 
der derivatives. These routines also work for f-functions. The original program's capabilities for 
the calculation of force-constants by numerical differentiation of gradients and for the optimiza¬ 
tion of geometries remain essentially intact. However there is now a choice of two optimization 
algorithms with the inclusion of Schlegel’s method. 

In addition to the above modifications a range of extra facilities were developed. These 
included a more powerful method of calculating one-electron properties and analyzing molecular 
charge distributions, and various ‘post-SCF’ stages, beginning with a 4-index transformation. 
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These new sections include M(f)ller-Plesset perturbation theory for total energies and molecular 
properties, and coupled Hartree-Fock calculations of molecular polarizabilities. The polarizability 
routines can calculate dynamic properties at real or imaginary frequencies, and obtain dispersion 
coefficients. It is also possible to use CHF theory to obtain the perturbations due to nuclear 
displacements. These can be used to obtain, analytically, all the dipole and quadrupole moment 
derivatives of a molecule. They also form part of the most important addition to the package, the 
section which calculates analytic second derivatives of the energy. This is a powerful technique 
whose speed and accuracy represents a considerable improvement over numerical differentiation. 

The capabilities of the latest CADPAC code (Version 4) include: 

a) the evaluation of one- and two-electron integrals over contracted cartesian gaussian 
basis functions of type s, p, d. or f. 

b) SCF calculations for closed shell, open-shell. UHF and generalized open-shell tech¬ 
niques. 

c) calculation of one-electron properties for these types of wavefunction, including a dis¬ 
tributed multipole analysis. 

d) calculation of the gradients of the SCF energy. 

e) use of the gradients for automatic geometry optimization, and for the calculation of 
force constants by numerical differentiation. There is a choice of two optimization al¬ 
gorithms. 

f) transformation of the integrals from the atomic orbital to the molecular orbital basis. 

g) M4>ller-Plesset perturbation theory calculations to third order in the energy and sec¬ 
ond order in the one-electron properties. 

h) coupled Hartree-Fock calculations of polarizabilities, including frequency dependence, 
and magnetizabilities. 

i) coupled Hartree-Fock calculations of the perturbation due to nuclear displacements. 

j) calculation of the dipole and quadrupole moment derivatives. 

k) calculation of the second derivatives (force constants) of the energy by analytic meth¬ 
ods. 

l) analytic calculations of polarizability derivatives. 

m) calculation of infrared and Raman intensities, and the study of vibrational circular 
dichroism. 
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n) calculation of MP2 gradients, dipole moment derivatives, polarizabilitites and force 
constants using analytic algorithms. 

o) spin-projected UHF MP2 energies 

2.2.4 HONDO 

The HONDO program, (Reference 46) originally developed by Dupuis and King at NRCC 
has recently been refined and updated by M. Dupuis at IBM-Kingston (Reference 48). 

The following features are available in the present version (HONDO 7) of the program: 

a) Single configuration self-consistent-field wavefunctions (closed shell RHF, spin unre¬ 
stricted UHF, restricted open shell ROHF), generalized valence bond GVB and gen¬ 
eral multiconfiguration self-consistent-field MCSCF wavefunctions, and configuration 
interaction Cl wavefunctions can be calculated. 

b) The electron correlation correction to the energy of closed shell RHF wavefunctions 
can be calculated by means of M4>ller-Plesset (MP) perturbation theory applied to 
second-, third-, and fourth-order (with or without the effects of triple excitations). 

c) The effective core potential approximation can be used. 

d) Optimization of molecular geometries using the gradient of the energy with respect to 
nuclear coordinates is possible with all but the Cl and MP wavefunctions. Optimiza¬ 
tion can be carried out in the cartesian space or in the internal coordinate space with 
the possibility of freezing some cartesian or internal coordinates. 

e) The force constant matrix in the cartesian space, and the vibrational spectrum, includ¬ 
ing infrared and Raman intensities, can be calculated with all but the Cl and MP 
wavefunctions. 

f) Calculation of the dipole moment and polarizability derivatives with respect to the nu¬ 
clear coordinates is possible, for use with a previously calculated force constant ma¬ 
trix. The force constant matrix can be transformed to the internal coordinate basis. 

g) Transition state structures can be determined with all but the Cl and MP wavefunc¬ 
tions by taking advantage of the energy gradients. 

h) The ‘Intrinsic Reaction Coordinate’ (IRC) pathway can be determined with all but the 
Cl and MP wavefunctions by taking advantage of the energy gradients. 

i) Molecular energies for several points on a potential energy surface can be calculated in 
a single run. 
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j) Non-gradient optimization of basis function exponents is possible. The source code 
can be modified to carry out optimization of other non-linear parameters, for example, 
contraction coefficients and geometrical parameters. 

k) The following electronic properties can be extracted from the wavefunction: 

1. dipole moment 

2. quadrupole moment 

3. Mulliken population, bond order and valency analyses 

4. spin density maps 

5. electron density maps 

6. electrostatic potential maps 

7. localized orbitals via Boys’ method 

8. static dipole polarizability 

9. static first and second hyperpolarizabilities. 

l) The potential due to finite point charges for a classical representation of an environ¬ 
ment, or a uniform electric field can be incorporated into the one-electron Hamil¬ 
tonian. 

2.2.5 COLUMBUS 

The COLUMBUS code (Reference 49 and 50) is a continuing development based on a joint 
project at Ohio State University (I. Shavitt), University of Karlsruhe (R. Ahlrichs), and Argonne 
National Laboratory (R. Shepard). The unique features of this code are the incorporation of 
relativistic or non-relativistic core potentials to permit analysis of heavy atom molecular systems 
and the inclusion of a sophisticated Cl package based on the unitary group approach of Shavitt 
(Reference 51). The basic programs included in the COLUMBUS code are as follows. 

1. AO Integrals 

This is R. Pitzer’s version of the integral package from HONDO with partial vectorization 
of auxiliary integral routines. This package can handle up to g-functions. 

2. Integral Transformation 

The transformation algorithm is written over all the orbitals which constitute a shell, using 
contraction coefficients defined over the primitive basis functions. Limited vectorization is 
possible over the innermost loops for dyadic operations. 
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3. SCF Gradients 

These routines are based on the HONDO version but include checks on integral symmetry 
to avoid operations over zero or near zero elements. At this time, only a first derivative 
analysis has been included. 

4. SCF Energy 

The SCF routines are typical of those found in the molecular structure codes. They incor¬ 
porate level shifting, damping and the incorporation of Pulay’s DfIS convergence accelera¬ 
tion procedure. 

5. MCSCF Analysis 

The MCSCF package is due to R. Shepard and incorporates extensive vectorization. The 
output vectors can be taken directly from this package and transferred to the Cl program. 

6. Multi reference Direct Cl 

This set of routines is based on the graphical representation of the unitary group ap¬ 
proach (GUGA) for constructing a Cl wavefunction. This code differs from conventional 
Cl in that a matrix representation of the hamiltonian is not explicitly computed. The uni¬ 
tary group approach can be much more efficient in most applications since ~80% of the 
hamiltonian matrix elements are usually zero but must still be included. This leads to the 
storage and diagonalization problems associated with large, sparse matrices. However, lit¬ 
tle work has been done on analysis of excited states with the same molecular symmetry as 
the ground state. An approach that looks promising is to define an approximate vector 
corresponding to the excited state wavefunction and to iteratively improve upon this solu¬ 
tion using the direct Cl contraction. The uncertainties in the procedure are applications to 
situations where degeneracies or near-degeneracies arise. Further studies of this case are 
in progress. 

2.3 Spin-Projected Unrestricted Hartree-Fock Method 

The unrestricted Hartree-Fock (UHF) method developed by Pople and Nesbet (Refer¬ 
ences 52 and 53) yields the best single-determinant approximation to the exact wavefunction for 
an atomic or molecular system. Such a wavefunction incorporates correlation by allowing orbitals 
of different spin to adjust to spatially different forms, thus breaking the symmetry restrictions of 
the conventional (RHF) method (References 54 and 55). It is necessary, however, to project from 
such a wavefunction, a properly antisymmetrized spin and angular momentum state in order to 
define eigenstates and eigenenergies corresponding to observable spectroscopic states. 

Let |a), refer to the doubly-occupied molecular orbitals (MO’s), (p a , <pp. Let fy), |<5), |r) 
refer to the singly occupied MO’s. We assume that all MO’s have been subjected to a 
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transformation to orthogonal form. Let be the coefficient associated with the permutation 
y «—► d (this permutation is (-1) times the overlap of the permuted spin eigenfunction with the 
original spin eigenfunction), and adopt the convention Cyy = -1. Then the expectation value of 
the nonrelativistic Hamiltonian 


K- 2 u <i) + XV(i,j) 

> <<j 


( 22 ) 


is given by 


(*) = 2X(«|U|a) + 2W U I)') + ZZM'W) 

a y aft 

- + 2 - XHvW (23) 

aft ay ay 

yi> yi 


where i, j refer to electron numbers, a, y, <5 etc., refer to MO’s and the sums run over all 

a, £.... 

The equations for determining the optimum MO’s can be derived from the following 
variational form 


6 (36)-2 {a\p)-Jj€ya{a\y) + {y\a))-'^€ yd (y|<5)l = 0 

. aft aft yd J 


(24) 
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These equations are 


(U + 2J C - Kc + Jo- —Ko)|a } = + -Z^£6a\d) 

fi 6 


(25) 


(U + 2JC-IQ + J 0 )|y)+ XC*l(<5| v |y)5) = ^ €flr \fi) + ) (26) 

<5 fid 


where 


|( d\V\y)d) » j dF0;(r)V(r,r)0.Xr)0 ( 5(r) 
(W|V|W') = f dFW + (r)V(r, P)W'(F) 


(27) 

Jc= 5>|V|0 Jo = X(<5|V|<5> 

fi d 

Kc|W = Xl(^l V l W )^) Ko|W = X!(<5|V|W)<5 ) 

fi d 


Left-multiplying Equation (25) by ( y| and Equation (26) by ( a| , these equations can be com¬ 
bined to yield 


^ = I( 2 C^+ 1)(*/<5|V|<5«) (28 ) 

d 


Noting that 


X Cyr(^ r IV | r/) = X Cdf(<5r | V | ry). 
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Equations (25) and (26) can be recast as follows. From Equation (25), 


L 6 

with 


(30) 


M = X( 2C <* + 1)|(3 }(r(<3|V|r)| 

dr 


From Equation (26), 


= X(2C /r + l)|^)(^r|V|(3 y) = O c M *\y) 
P P x 


with 


o,-£l/J)(/»|; Oo - Z I* 5 K**! 

P <> 


(32) 


Using Equations (31) and (32) in (25) and (26), we find 


(U + 2J c -Kc +J 0 -|Ko-^M)|a)= £(f,a\P) 


(33) 


(U + 2J c -Kc + Jo-^Ko + ^M t -O c M t )|y)= £e*|<3 
1 6 


(34) 
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Finally, Equations (33) and (34) can be rearranged to yield the same Hermitian operator on the 
left-hand side with the resulting definition of the one-electron eigenenergies for the closed and 
open-shell eigenstates: 

€a' =(a|U + Uc-fQ + J 0 -^Ko|a ) (35) 

e/ =(y|U + 2J c -Kc + J 0 -jK 0 + |M t |y) (36) 


Equations (35) and (36) represent the correct one-electron energy expressions for the spin- 
projected eigenstate. 

This formalism has been developed to extend the conventional HF method to include 
split-shell correlation and proper spin and symmetry projections. It is being incorporated into a 
computer program using gaussian-type orbitals (GTO’s) as the eltmentary basis functions. A 
crucial feature of this method is that dissociation always follows the lowest energy pathway 
thereby permitting a proper description of bond formation and breakage. In contrast. RHF or 
MC-SCF methods often exhibit improper dissociation character or exhibit size inconsistencies 
owing to correlation energy changes in going from molecular geometries to separated atom- 
molecule or atom-atom dimensions. 

2.4 Transition Probabilities 

The electronic and vibrational-rotational wavefunctions of a pair of states can be used to 
calculate transition probabilities. If two molecular states are separated in energy by an amount 
AE nm = hcv (h = Planck’s constant, c = velocity of light, v = frequency in wave numbers), the 
semi-classical theory of radiation (References 57 and 58) yields for the probability of a spontane¬ 
ous transition from an upper state n to a lower state m 


Ann, « 


4 AEL S, 


nm J nm 


3 tC C 3 g n 


(37) 


Here A nm is the Einstein coefficient for spontaneous transition from level n -+ m, g n is the total 
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degeneracy factor for the upper state 


g n = (2 - d 0 , A')(2S’ + 1)(2J' + 1) 


(38) 


and S nm is the total strength of a component line in a specific state of polarization and propa¬ 
gated in a fixed direction. A related quantity is the mean radiative lifetime of state n defined by 




m < n 


(39) 


the summation being over all lower levels which offer allowed compactions. The intensity of the 
emitted radiation is 


Inm — ^E nm N n A n 


(40) 


where N n is the number density in the upper state n. This analysis assumes that all degenerate 
states at the same level n are equally populated, which will be true for isotropic excitation. The 
total line strength S nm can be written as the square of the transition moment summed over all de¬ 
generate components of the molecular states n and m: 


S 


nm “ 



(41) 


where j and i refer to all quantum numbers associated collectively with upper and lower electronic 
states, respectively. 

In the Bom-Oppenheimer approximation, assuming the separability of electronic and 
nuclear motion, the wavefunction for a diatomic molecule can be written as 


^'vjma = ^’ci(r.R)vvfR)V9MA(0.;r.0) 


(42) 
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where ^(F, R) is an electronic wavefunction for state i at fixed internuclear separation R, tp v { R) 
is a vibrational wavefunction for level v, and i/jjam (#,*.0) refers to the rotational state specified 

by electronic angular momentum A , total angular momentum J and magnetic quantum number 
M. The representation is in a coordinate system related to a space-fixed system by the Eulerian 
angles ( 9,x,<p ) • The transition moment Mj, can be written, using the wavefunction given by Equa¬ 
tion (42), as 


M ji = J Vvj'A'M'{^ e + "A’ M ^e^TvJrr (43) 

The subscripts e, v and r refer to the electronic, vibrational and rotational wavefunctions and M e 

and M n are the electronic and nuclear electric dipole moments, respectively. Integration over the 
electronic wavefunction. in the Born-Oppenheimei approximation, causes the contribution of the 

nuclear moment M n to vanish for i # j. The electronic dipole moment can be written (Refer¬ 
ences 57 and 58) in the form 


= -Z eFk ' = ~{ X e? k} (44) 

k k 

where the primed coordinates refer to the space fixed system, the coordinates Fk refer to a mole¬ 
cule-fixed system and ®{0,x, (p) is a group rotation tensor whose elements are the direction co¬ 
sines related to the Eulerian rotation angles (9rx, <p >). Using bracket notation. Equations (43) and 
(44) can be combined to yield for the transition moment 

Mj» = M 1 v' J Av'M" =( jM -X ef k|iv") •( J'A'M’|i(0,* ( 0)|J"A”M'') ( 45 ) 

k 

The matrix elements ( J'A'M’ |‘35(0,^,0)|J"A"M”) determine the group selection rules for an 
allowed transition and have been evaluated for many types of transitions (References 59-61). Sum¬ 
ming Equation (45) over the degenerate magnetic quantum numbers M' and M” we have from 
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Equation (41) 


q _ Qnv'JW _ *J'A' ^nv’ (^\ 

3 nm ~ *>mv"J"A" “ J J"A”Pmv'' 

where is the Honl - London factor References (62 and 63) and 

Pmv • = X K J v ' I - X ef k I iv ”) I' 2 (47) 

i.j k 

is the band strength for the transition. Combining Equations (38), (40) and (46), we have for the 
intensity of a single emitting line from upper level n: 


I 


nm — 


;nv'J 

'mv’T 1 


, fACnv'J' l^cnv'J'A' 

4 [ ac miVJ”J a mv "J"A' 

3 J ' /j 4 c 3 w n (2f' + 1) 


(48) 


where Nj' is the number density in the upper rotational state J' and to n = (2-<5 oA ) (2S' + 1) is 
the electronic degeneracy. Taking an average value of AE^ V V »V- for the whole band. Equation (48) 
can be summed to yield the total intensity in the (v', v') band: 


I 


nv' 

mv" 


- I CV 

I' I" 


a TACnv' l 4 r,nv' 

4 , Pmv' - 

T >4 3 

3 h e (On 


(49) 


where N v ' = ^ Nj> is the total number density in the upper vibrational level v' and where we 
y 

make use of the group summation property 


I*rX" = (21' + 1) (50) 

I" 


Comparing Equations (40) and (49) we have for the Einstein spontaneous transition coefficient of 
the band (v', v") 


A 


nv' 

mv" 


4 [AE£..fp£.. 
3 h A c y (o„ 


(51) 
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Similarly, the lifetime of an upper vibrational level v' of state n can be written 


7- XI A- ,52) 

n m<n v" 

where the summation runs over all v" for each lower state m. Equation (51) can be cast in the 
com^utati^nui form 




(21.41759xl0 9 ) 


AE! 


U) n 


“ V ..(a.u.)] 3 pmv(a.u.l 


(53) 


where AEJ}, V V .. and p^ v v -. are in atomic units. It is also often convenient to relate the transition 
probability to the number of dispersion electrons needed to explain the emission strength classi¬ 
cally. This number, the f-number or oscillator strength for emission, is given by 


mc 3 h 2 


2e 2 [AE™ v .; 


(54) 


The inverse process of absorption is related to the above development through the Einstein 
B coefficient. Corresponding to Equation (40), we have for a single line in absorption 

= | K(v)dv = hv mn N m B mn 

line(v"v’J”J') 

where K(v) is the absorption coefficient of a beam of photons of frequency v and 


B 


mn 


= B 


nv’J'A' 

mV'J"A" 


2ji 


cnv'J'A' 

^mv'T'A' 


3fl 2 C 0) m [ 2J" + 1) 


(56) 


is the Einstein absorption coefficient for a single line. Summing over all lines in the band (v ", v ), 
assuming an average band frequency, we obtain 


= N v 


271 n nv '..AF n 
rmv 


jh 2 CO) n 


(57) 
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where N v - = ^Nj- is the total number density in the lower vibrational state v". Corresponding 
r 

to Equations (53) and (54) we can define an f-number or oscillator strength for absorption as 


f - p nv 

In computational form. Equation (58) becomes 

f . 2 . AEt-(a-U-) y 

‘mn.v' v - - p mv -ta.UJ (59 

j w m 

where AE™... and pj^ v .. are in atomic units. Combining Equations (51) and (54) and comparing 
with Equation (58), we see that the absorption and emission f-numbers are related by 


l mn,v v — V 7*nm,vV 
(O m 


Some caution must be observed in the use of f-numbers given either by Equations (54) or (58, 
since both band f-numbers and system f-numbers are defined in the literature. The confusion 
arises from the several possible band averaging schemes that can be identified. 

An integrated absorption coefficient (density corrected) can be defined from Equation (57) 


Sv-y = A| £ .. = N V B, .(|-exp ~ V )h~^p I 

where the exponential factor corrects for stimulated emission. Equation (61) can be written in 
terms of the absorption f-number as 


7re 2 Nv''/, -hcvv'-.v'x, 

Sv " v ' “ me 2 P ( CXP kT )fmn v 


32 






Using hc/k = 1.43880 cm-K. we obtain a computational formula for the integrated absorption 
coefficient as Sy • V '(cm" 2 atm' 1 ) = 


23795 * 107 <S^<l-exp 14388 °--^ - > ) U . 


i(K°) a N t 

The total integrated absorption is found from 


(63) 


Stotal XX^ v ”- v ' (64) 

v" v' 

where, unaer normal temperature conditions, only the first few fundamentals and overtones con¬ 
tribute to the summations. 

The developments given above are rigorous for band systems where an average band 
frequency can be meaningfully defined. Further approximations, however, are often made. For 
example, the electronic component of the dipole transition moment can be defined as 


%(R) =(jl -2er k |i) (65) 

k 

This quantity is often a slowly varying function of R and an average value can sometimes be cho¬ 
sen. Equation (47) can then be written approximately in factored form as 

Pmv" “ fivV' (66) 

•0 

where q V V ", the square of the vibrational overlap integral, is called the Franck-Condon factor. 

% is evaluated at some mean value of the internuclear separation R. In addition, it is sometimes 

possible to account for a weak R-dependence in M e by a Taylor series expansion of this quantity 
about some reference value, R^ usually referred to the (0, 0) band. We have, 

% => %jf[l + a(R-R^) + b(R-R^) 2 + ...] (67) 
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Substituting into Equation (67) and integrating yields 


Pmv" - qw-X l<* [ l + a(Rw-R^) + b(R vV - R,*) 2 + ...] | 2 


n >—ft - \ v 'i( R -Miv") 

(K v v'-K^-7-777771- 


is the R-centroid for the transition and 

' (V | v ') 

is the R 2 -centroid. Note that this last term differs (to second order) from the square of the R-cen¬ 
troid. An alternate procedure can be developed by evaluating Equation (65) at each R-centroid, 
R v v ■ Then 

Pmv" ~ q»V" 1) 

Equation (71) assumes that the vibrational wavefunction product tp v - . behaves like a delta 
function upon integration, 

WtfV- = d(R-R‘vv)(v'|v'*) (72) 


The range of validity of Equation (71) is therefore questionable, particularly for band systems 
with bad overlap conditions such as oxygen Schumann-Runge. The range of validity of the R-cen- 
troid approximation has been examined by Frazer (Reference 64). 

The final step in calculating transition probabilities is the determination of $>,,(R), the 
electronic dipole transition moment, for the entire range of internuclear separations, R, reached in 
the vibrational levels to be considered. This can be expressed in terms of a Cl expansion as 

&).( R ) = ) (73) 

liV 

where ci and ci, are coefficients for i// eI and ij/ c[ , respectively. 
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An analysis similar to that yielding Equations (10) and (11) gives 


n n 

(^(R)|M c |^v(R)} = 2 e p(^slOsP|^M s )( n^fr-R)|M e P| n Mnc.R)) (74) 

p k = 1 k = 1 

The spatial integral in Equation (74) reduces to one-electron integrals equivalent to overlap inte¬ 
grals, and the evaluation of Equation (74) can be carried out with any of the standard ub initio 
programs such as DIATOM or GAMESS. Programs for evaluating 3bj,(R) in Equation (73) have 

been developed at UTRC and examples of their application have appeared in the literature (Ref¬ 
erence 65). 

For perturbed electronic systems, the transition dipole moment will have a strong R-depend- 
ence and R-centroid or other approximations will be invalid. A direct evaluation of Equation (47) 
would therefore be required using the fully-coupled system of electronic and vibrational 
wavefunctions to properly account for the source of the band perturbations. 
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SECTION 3 


DISCUSSION AND RECOMMENDATIONS 

Several broad surveys of advanced propulsion systems have been undertaken within the past 
ten years (References 66-70). The objectives of these surveys were to identify new sources of 
propulsion energy, ranging from novel energetic chemical species to advanced concepts such as 
antimatter annihilation (References 71 and 72). Within the framework of advanced propulsion 
systems deriving their energy source from non-nuclear chemical reactions, we can consider several 
categories of energetic species. They include: 

1 ) ordered (or metallic) forms of matter which at conventional temperatures and pres¬ 
sures are molecular solids. Metallic hydrogen and ammonium are examples. 

2 ) free radicals, either atomic or molecular in nature, that are stabilized against recombi¬ 
nation by condensation and immobilization on noble gas matrices held at very low 

(4 K) temperatures. Atomic hydrogen, beryllium monohydride (BeH) and the hydroxyl 
(OH) radical are examples in this categoiy. 

3) electronically excited metastable atoms or molecules that are stabilized against radia¬ 
tive decay by externally applied electric, magnetic or radiation fields, or by an intrinsic 
field within a condensed phase of such material which acts as a barrier to their decay. 
Metastable atoms, such as He"( 3 S) or N*( 2 D), and metastable molecules, such as 
Hjf'Ai), He 2 (a 3 Zu or 5 ), HeH^A'II) and H^B 1 !^), are examples in this cate¬ 
gory that would constitute significant energy sources (> 1000 kjoule/mol). 

4) strained molecules with high positive heats of formation. This is an important class of 
molecules which are really isomeric forms of conventional chemical species but are 
structurally prevented from relaxation to their lowest energy state or to decomposition 
by significant activation barriers on their potential energy surface. Examples in this 
category include vinylidine, several structures of Li m H„ , cyclopropene, aziridine, 
azetidine and covalent azides or azide-like structures. Hypervalent compounds such as 
NF 5 or CIOF 5 would also be included in this category. 

Ordered or metallic states of solids that are normally molecular have been extensively 
studied both experimentally and theoretically (References 73-76). Detailed quantum mechanical 
studies of several crystalline forms of metallic hydrogen have been carried out which yield 200-300 
kjoules of energy above the normal solid molecular form. This translates into an available specific 
impulse in excess of 1200 seconds. However, the derived thermodynamic properties from these 
calculations suggest that external pressures of the order of several megabars would be required to 
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form the metallic phase. Since this pressure requirement is approximately one order of magnitude 
higher than has ever been achieved (except possibly by thermonuclear explosions), there is 
currently very little active research in this area. There appears to be little more that theoretical 
studies can contribute until significant advances in high pressure technology are achieved. 

The concept of using stabilized free radicals as a propulsion energy source must be tem¬ 
pered by an assessment of the maximum specific impulse that can be derived from such systems. 
The bond energy of molecular hydrogen (436 kjoule) translates into an available specific impulse 
in excess of 1800 seconds, provided there are no constraints on the density of hydrogen atoms 
that could be formed and stabilized. The available specific impulse drops dramatically when other 
free radicals are considered. As an example, atomic oxygen atoms, stabilized in their ground 3 P 
state, would yield approximately 600 seconds and most other dissociated molecules, stabilized as 
free radicals, are not sufficiently energetic to yield a specific impulse greater than the baseline 
value of 389 seconds available from H 2 /O 2 . The early research efforts in free radical stabilization 
are reviewed in the book by Bass and Broida (Reference 10). An analysis by Rosen (Reference 77) 
on the prospects for stabilizing hydrogen atoms in significant concentrations, and the recent 
experimental studies on spin polarized hydrogen reported by Hess (Reference 78), suggest that 
storage of high densities of hydrogen atoms will be difficult to achieve. The very low temperatures 
required (< 4 K) to minimize recombination processes and the dilute surface storage require¬ 
ments are serious obstacles when viewing these materials as chemical propellants. An experimen¬ 
tal program was carried out for AFRPL by this Center in an effort to stabilize light element 
atoms at 77 K (Reference 79). The maximum concentration that was realized was well below that 
required for consideration as a possible propulsion system. Attempts to stabilize H atoms in 
various host matrices all indicated that no significant reaction barrier to atomic recombination 
could be achieved. These observations, coupled with a thermodynamic analysis indicating that 
only stabilized light element atoms would offer any significant improvement over conventional fuel 
types, suggest that this area may be of limited value for consideration as an advanced fuel type. 
Recently, however, (References 80-82), the idea of trapping ions or free radicals in a noble gas 
matrix has been reexamined under the High Energy Density Matter (HEDM) program. The basic 
concept is that noble gas hydride species such as Ar„H + are stable and may offer a route to 
hydrogen atom stabilization. Another possibility is to store metal atoms interstitially in a dilute 
host matrix such as hydrogen. The intermolecular separation of H 2 molecules in the solid phase 
is large compared to many light metal atom radii. Both theoretical and experimental studies of 
these mechanisms should be further pursued. 

The above analysis, however, suggests that research efforts for identifying high energy 
non-nuclear fuels should probably be concentrated in the remaining categories of metastable 
species and highly strained molecular compounds. Since the number and type of compounds that 


37 



might be formed is potentially very large, it would be helpful if certain criteria could be estab¬ 
lished to focus both theoretical and experimental studies into areas where there is a significant 
higher probability of forming high energy molecular species. A systematic study to establish such 
criteria for screening diatomic and triatomic molecules that could be formed from the first 18 
elements was performed by this Center for AFRPL in the mid-sixties (Reference 83). The 
electronic structure and energy content of those molecules that showed high performance 
potential as energy sources in a chemical propulsion system were then examined using state-of- 
the-art quantum mechanical valence bond programs. These calculations were pioneering at the 
time but are admittedly crude in comparison to results obtainable today from currently available 
techniques, as outlined in the Mathematical Background section of this report. 

This research program for AFRPL, which surveyed some 2000 electronic states of 171 
diatomic molecules and over 10,000 electronic states of 324 triatomic molecules that could be 
formed using atoms ranging from H to Ar, did establish certain fundamental criteria for estimat¬ 
ing performance potential. The foremost guideline involves the molecular weight of the energetic 
species. The specific impulse of chemical propulsion systems can be estimated on the basis of the 
properties of the major exhaust products, i.e. the effective species which determine the available 
energy content of the reaction. To a good approximation we Find: 



where I sp is the specific impulse (sec), AH is the difference between the combustion chamber and 
exhaust enthalpies (kcal) and W is the weight of the propellant (gm), as determined by the effec¬ 
tive number of moles of gaseous products and the stoichiometry of the propellant system. For 
many propellant combinations, it has been shown that the quantity, AH/n e , which represents the 
change of enthalpy per mole of gas phase species, can be related to the heat of reaction of the 
propellant or, more simply, to the available energy content. If hydrogen gas is chosen as the work¬ 
ing fluid, a universal plot of specific impulse as a function of molecular weight of the energetic 
species, with energy content as a parameter, can be constructed. This is illustrated in Figure 1 
where we have also labeled the specific impulse available from both conventional fuels and those 
proposed as high energy fuels. 

It is clear from Figure 1 that the molecular weight of the ‘‘energy carrying” species in a 
chemical propulsion system is the overriding factor for our evaluation of expected performance. 
Taking metastable He 2 (a 3 Z^ ] as an example (labeled 15 on Figure 1), we find a maximum 
available I 5p of ~ 1370 seconds. The chemically similar metastable NeJ, however, only yields a 
maximum I, p of ~ 630 seconds, mainly due to the detrimental effects of increased molecular 
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SPECIFIC IMPULSE, Isp (sec) 


HYDROGEN CONTENT VARIED TO OPTIMIZE IMPULSE (1000-14.7 psia) 



MOLECULAR WEIGHT OF ENERGETIC SPECIES 


Figure 1. Specific Impulse Obtainable from Fuel/Oxidant Energy Release 
with Hydrogen as a Working Fluid 
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weight. Those species labeled as high energy fuels on Figure 1 are all capable of delivering specific 
impulses greater than ~ 1000 seconds except for Ne 2 and U 3 H*, which is just at baseline I sp 
(450 sec), again mainly from molecular weight considerations. Using an available energy content 
of ~ 1200 kjoules which, except for Hej, is a realistic upper limit to the energy that can be stored 
in an electronically excited metastable state of a molecule, we find a molecular weight cut-off of 
~ 30 if a 50% improvement in performance is desired over that available from conventional fuel 
types. This is obviously not a hard number for all metastable species but gives us guidance to the 
type of compounds that should be considered. 

In the category of strained molecules with positive heats of formation, we can establish an 
appropriate upper limit to the available energy of such compounds based on their lowest 
ionization potential. The reasoning here is that isomeric forms of a given compound, such as 
U 3 H, which have energy in excess of their lowest ionization potential, would constitute autoioniz- 
ing or continuum electronic states with lifetimes dictated by intrinsic quantum mechanical decay 
pathways (spin-recoupling, electron autodetachment, etc.). Considering only compounds formed 
from atoms through carbon, we again find ( He being an exception) a molecular weight cut-off of 
~ 30 for a 50% improvement in I sp over conventional fuels. A more conservative (and realistic) 
estimate of the molecular ionization potentials (600 kjoules) would indicate that emphasis should 
be placed on studies of energetic species with molecular weight of 20 or less. This clearly rules out 
all hydrocarbons and even most boron compounds if significant increases in chemical propulsion 
performance are desired. We are thus led to consider, as primary candidates , those molecular 
structures that can be formed from hydrogen through boron which have a molecular weight of 
~ 30 or less. The chemistry of boron compounds has been extensively studied (References 84 and 
85) and only those energetic compounds containing at most three boron atoms would meet our 
criteria for a significant improvement in available specific impulse. 

This preliminary analysis does not rule out the possible utility of higher molecular weight 
species as advanced fuel types. Our studies of both advanced fuels and oxidizers, under this 
research program, have indicated (to date) that the identification of an advanced fuel with 
performance significantly better than H 2 • using O 2 as an oxidizer, may prove to be difficult. The 
C 2 V structure of, U 3 H may offer some improvement and is described in more detail below. 
However, our studies of oxidizers have resulted in the identification of the covalent azide. FN 3 , 
and the azide-like molecule, a-N 202 . as powerful new oxidizers when used as replacement 
species for ammonium perchlorate (AP) in typical solid propellant formulations. In addition, 
theoretical studies indicate that the hypervalent structure, NF 5 , is vibrationally stable with a 
positive heat of formation for the D 31 , (trigonal bipyramid) structure. Although previous attempts 
to synthesize this molecule have been unsuccessful, our calculations shed some light on synthesis 
routes that may be promising. Our preliminary studies of high energy species have thus led us to 
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the conclusion that significant improvements in solid formulations may be possible through the 
development of advanced oxidizers. 

Turning next to an analysis of electronically excited metastable molecules, we note that a 
theoretical study of such species has been performed by this Center for AFRPL (Reference 83). 
This study and more recent investigations acknowledge that no mechanism has yet been found 
that would prevent radiative decay of such molecules in a time span of a few seconds or less. 
Recently, however, Zmuidzinas has proposed (Reference 13) that spin-polarized triplet helium 

atoms could be stabilized by forming the metastabie 5 I g + state of Hei which exhibits no first 

order decay process to the autoionizing continuum of He^ + e or to ground state helium atoms. 

In an effort to more fully understand the chemistry of such excited states we lccentlv 
initiated a series of calculations of the low-lying quintet states of He 2 using extended Cl 
wavefunctions. A previous ab initio study by Garrison, Miller and Schaefer (Reference 15) had 
indicated that the molecular state of He 2 , arising from the interaction of two 3 S(2s)He* 

atoms, exhibits 0.07 eV stability with no long range repulsive barrier. In Figure 2 we outline the 
results of our recent study. The low-lying 5 Ig quintet state of He? arises from the 

He* 3 S(2s) + He* 3 S(2s) asymptote while the low-lying , 5 n u , 5 H g and 5 £+11 states 
arise from the He* 3 S(2s) + He* 3 P(2p) asymptote. All of these molecular states are modified by 
interactions with the quintet 5 Zg , 5 I u + . 5 n g and 5 n u components of the ionic states arising 

from the He'[ 4 P°(2s2p)] + He + [ 2 S(2s)] asymptote. We find that the 5 H K state is actually the 
lowest molecular state among this group and would act to completely predissociate the weakly 
bound 5 Ig state. A further analysis of possible decay mechanisms for this low-lying 5 n g state 
indicates that a maximum radiative lifetime of ~ 10 n sec could be expected for this metastable 
system. Thus stabilization or metastable He*[ 3 S(2s)] atoms by quintet molecular state formation 
is no more likely than by formation of the lowest a 3 S^ excited state of He 2 . 

We are led to the conclusion that simple stabilization of electronically excited atoms through 
molecule formation is unlikely. However, it might be possible to utilize the changed chemical 
characteristics of an excited atom or molecule to form a new molecular structure which represents 
a local point of stability on the ground state potential energy surface. Illustratively, we have the 
following example in Figure 3. 
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INTERNUCLEAR SEPARATION R (A) 


ENERGY SCALE IS GIVEN IN ELECTRON VOLTS (1 eV-96 54 k|Oules/mol) AND IS RELATIVE TO 

He [ 1 S (1S 2 )) + He (’S (is 2 )]-0. THE IONIC CURVE LABELED BY THE SQUARES HAS THE He” ( 4 P] + He* | 2 S] 

DISSOCIATION ASYMPTOTE. 


Figure 2. Low-Lying Quintet States of He 2 ' 
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Figure 3. Generalized Reaction Surface 

Here the most stable molecular configuration corresponds to £1 whereas the system at £2 repre¬ 
sents an isomeric form of the molecule with extra energy, E ex and which is stable against decay 
to the lowest state by the presence of an activation barrier, E a . The crucial questions are: how 
high is the activation barrier in any system and how much extra energy can be stored in the sys¬ 
tem stabilized at this local minimum. The loss of molecules from the system at £2 can be realized 
by tunneling (if E a is small) or by extra internal energy as vibrations or rotations which would 
exceed E a . The ideal system would have a large E ex (high energy storage) and a large E a (high 
barrier for decay) with one or more vibrational levels possible at the geometry represented by £2 ■ 
It is tempting to call these structures isomeric conformations of the most stable ground state ar¬ 
rangement of atoms, but that label is usually reserved for regions £2 that can be reached by sym¬ 
metry operations (rotations, reflections) of the molecule. Often the molecular conformation at 
£2 represents a high symmetry point on the potential energy surface. It is crucial to examine such 
points to determine if they represent true minima or if they are saddle points which are unstable 
against some mode of internal molecular motion. We show below that the H 4 (tetrahydrogen) sys¬ 
tem is a classic example of a local minima on a ground potential energy surface which, in fact, 
represents an unstable saddle region. 

Structural and thermodynamic studies of molecules which correspond to conformations 
represented by the region £ 2 have been carried out under the present contract. To date, calcula¬ 
tions have mainly focused on four types of compounds: 

1) light element Cj v and C 2 v structures (H 4 , U 3 H, LiHj, Li 4 ) which represent ad¬ 
vanced fuel types 

2) halogen azides and azide-like structures (FN 3 , FNCO, C0 3 , HN 3 , FNBF, a-N 2 0 2 ) 
which represent advanced oxidizers 
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3) cyclic boron structures (B 3 H 3 , B 2 H 2 NH) which are analogous to strained ring hydro¬ 
carbons cycleproplyene (C 3 H 4 )and ethyleneimine (C 2 H 5 N) 

4) hypervalent molecules (NF 5 , PF 5 ), which exhibit high oxidization states. 

We summarize below our theoretical studies in these four areas and suggest new avenues of 
research based on our experience under the present contract. 

3.1 Light Element High Symmetry Structures 

3.1.1 H 4 

Our initial theoretical studies of advanced fuel types were focused on the energy content and 
stability of the H 4 (tetrahydrogen) system. The simple H 4 system has been the subject of many 
quantum chemical studies (Reference 86 ). Until recently, it was thought that the only binding in 
the ground state was due to the weak van der Waals interaction of two H 2 molecules (Refer¬ 
ence 87) although several workers have suggested the existence of bound excited states (Refer¬ 
ence 88 ). In particular, Nicolaides, Theodorakopoulos, and Petsalakis (References 89 and 90) have 
predicted the existence of a 1 A' bound excited state of H 4 in C s symmetry, formed from the 
interaction of ground state H 2 with vibrationally excited He 2 B 1 ^ . The calculations of 
Nicolaides, et al. show that this bound excited state of H 4 arises from an avoided crossing with 
the lowest 1 A' state. It is therefore likely that this first excited 1 A' state of H 4 is predis¬ 
sociated through vibronic interaction with the ground state. In pure C 3v symmetry, the two 1 A' 
states discussed by Nicolaides, et al. correlate to the ground 1 Aj state and an excited 1 E state. 
Because these states are of different symmetry, there is no avoided crossing in C 3 V symmetry to 
stabilize the excited state, as there can be in the lower C s symmetry. However, these same 
calculations also indicated that the ground state surface exhibits a local minimum in C 3 V 
symmetry. Physically, this conformation may be viewed as arising from the binding of a hydride 
ion (H") to the stable trigonal ion. The purpose of our initial study was to examine this 
conformation in C 3 V symmetry to determine the character of this minimum point. 

Ab initio calculations on H 4 were performed with the 6slp/4slp basis of Nicolaides et al., 
and a more accurate 10s3p/5s3p basis. The I0s3p/5s3p basis consists of the 9s/4s contraction of 
Siegbahn and Liu (Reference 91) augmented by a diffuse s-type function (exponent 0.016). Using 
the Gaussian 82 system of programs (Reference 41), C 3 V optimized structures were found at the 
SCF, MP2, and CISD levels of theory using the 6slp/4slp basis, and at the SCF and MP2 levels 
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using the 10s3p/5s3p basis. At each optimized geometry, a vibrational frequency analysis was 
performed. A summary of the results of these calculations is given in Table 1. 

The most significant finding of these calculations was the existence of two imaginary 
vibrational frequencies, corresponding to an unstable vibrational mode of e symmetry. This result 
was obtained with both large and small basis sets, and at every level of theory. Using the 
6slp/4slp basis set, an SCF calculation of the path of steepest descent from the C 3 V optimized 
geometry was performed with the intrinsic reaction coordinate pathfinder of Schmidt, Gordon, 
and Dupuis (Reference 92) implemented in the GAMESS code. This calculation was performed 
with no symmetry constraints and indicated a monotonically decreasing path which dissociates 
into two H 2 molecules. We have concluded that the H 4 energy minima found in C3,, symmetry 
are not stable with respect to geometric perturbations which lower the spatial symmetry and that 
the ground 1 Ai state of H 4 in Cj v symmetry is not bound. Exploratory calculations were also 
performed in C s symmetry to ascertain the possible existence of a bound state of H 4 in a slightly 
distorted pyramidal configuration. We were unable to locate such a state, as all of our optimiza¬ 
tions tended either to dissociation or to the C 3 V minimum described above. 

Our stability analysis of H 4 (C 3 V ) was obviously a disappointing result since, as shown on 
Figure 1, such a molecular species is estimated to yield a specific impulse of nearly 1300 seconds, 
based on our calculated bond energies relative to two ground state hydrogen molecules. Other 
structures of this type, however, are suggested by the similarity of H~ and Li' . In Figure 4 we 
indicate several potential high-energy, light element C 3 V molecular structures that result from 
systematic substitution of Li for H atoms. Li' has an electron affinity of 0.62 eV. very similar to 
that of H" (0.75 eV). In addition, Li 3 forms a stable molecule (Reference 93) with a binding 
energy of ~ 1.5 eV relative to ground state Li + Li ’2 . whereas the neutral H 3 system is unstable 
at all geometries. This suggests a greater stability for the Li 3 H species, but as shown in Figure 1, 
the increased molecular weight causes a penalty in the predicted I, p . 

3.1.2 Li 3 H 

Based on this analysis, the C 3 V structure of L 13 H, which is the analog of tetrahydrogen with 
a Li^ base, was examined using several basis sets and at the SCF and MP2 levels of theory. 
Calculations at both the SCF and MP2 level of theory, using several basis sets including diffuse 
orbitals are given in Table 2. These data indicate that this species is stable in all approximations 
that were examined. A diffuse basis (6-31 ++ G *) was examined at the SCF level to insure 
that Li* and H' were accurately represented. However, the diffuse functions appear to have little 
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Table 1 . Ab Initio Structure and Calculated Energy for H4 in C 3v Symmetry 


Theory 

h 

r 

Energy 


Frequencies 






e 

a 1 

e 

a\ 

SCF/small 

1.8784 

0.8433 

-1.982 831 

2942/ 

1053 

3069 

3794 

MP2/small 

1.8094 

0.8477 

-2.035 124 

3614/ 

1096 

3010 

3751 

CISD/small 

1.8048 

0.8509 

-2.049 161 

4025/ 

1092 

2962 

3712 

SCF/large 

1.8634 

0.8421 

-2.000 761 

2188/' 

1167 

3101 

3777 

MP2/large 

1.8170 

0.8451 

-2.063 595 

2503/ 

1201 

3051 

3740 


Energies (hartrees), distances (A), frequencies (cm -1 ) 

h is the distance of the fourth proton from the center of the H 3 triangle of side r. 

Small and large refer to the 6slp/4slp and 10s3p/5s3p basis sets, respectively. Note that the fre¬ 
quencies of the e modes are doubly degenerate. 
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c) LigH (STABLE) 


d) Li 4 (UNSTABLE) 


Figure 4. High-Energy Light Element C 3 V Molecular Structures 







Table 2 . Theoretical Results for U3H C3 V Structure 
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effect in geometry optimization. The results for this basis were essentially the same as those with 
the double-zeta plus polarization basis sets. The optimized structure for U 3 H has a short 

pyramid height (~ 1.1 A ) for the H' anion above the U 3 base. The structure thus somewhat 
resembles the NH 3 molecule. Although it is tempting to assign the ion pair state (Li^ + H‘) 
interaction to the trigonal pyramidal structure (C 3 ,,), because of its structural resemblance to the 
carefully studied H 4 species, the disparity in the Li—Li distance (at MP2/6-311 ++ G * *) be¬ 
tween this structure (2.577A) and Li^ (3.001 A) suggests that the C 3 V structure is not correctly 
characterized as an Li^ + H" ion pair. This disparity was also noted by Talbi and Saxon 
(Reference 94), who suggested that the pyramidal structure is of mixed ion-pair and covalent 
character based on a wavefunction analysis. 

Recently Cardelino, et al. (Reference 95) have shown that a kite-like C 2 V structure for U 3 H 
is the ground state at the SCF level of theory. Calculations of this structure were carried out at 
both the SCF and MP2 level of theory, using 6-31G* and 6-311G** basis sets. The data are 
shown in Table 3 where the C 2 v planar kite-like structure is, in fact, found to be lower in energy 
than the trigonal C 3 V structure at both the SCF and MP2 level of theory. A harmonic frequency 
analysis indicates vibrational stability for this C^ structure, which appears to be a Li atom 
bridged across the known, stable L^H molecule. The ionicity on the H atom is ~ 18% in this 
conformation as compared with ~ 29% for the trigonal C 3 V structure. Since both of these 
structures exhibit stable vibrational frequencies, they probably represent different local minima on 
the same ground state surface. Preliminary studies of the transition state connecting these two 
structures have been carried out which indicate a slightly distorted (C s ) trigonal bipyramid, with 
nearly the same energy as the C 3 V structure. Further detailed calculations of L^H are given in 
Appendix C where we illustrate that inclusion of zero-point energy requirements places the C 3 V 
structure at the same level as the transition state, making the thermodynamic stability of the C 3 V 
structure questionable. However, the C*, structure is found to be stable, and represents the 
absolute ground state minima of the L^H surface. 

The performance of the Li 3 H/F 2 and IJ 3 H/O 2 bipropellant systems has been evaluated 
through specific impulse calculations. Calculations were performed for a combustion pressure of 
800 psia and equilibrium flow through a nozzle having an exit-to-throat area ratio of 20:1. Two 
sets of calculations were performed in which the L^H was taken as the C 2 v structure 
(AHf = -3.49 kcal/mol) and as the C 3 V structure (AHf = + 23.83 kcal/mol). As noted above. 
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Table 3 . Theoretical Results for Li 3 H C 2v Structure 
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only the C 2 V structure appears to be vibrationally stable. The calculated performance curves are 
given in Figures 5 and 6 . respectively, using F 2 and O 2 as oxidizers. It is seen that bipropellant 
combinations offer approximately a 10% improvement over the UH/O 2 or LiH/F 2 systems. 
However, a significant increase in performance is available when U 3 H is used in a tripropellant 
formulation. The peak I sp using O 2 and an optimum H 2 ratio (IJ 3 H/H 2 = 0.5), and assuming 
the C 2 v structure for L 13 H, is 456 sec., which represents a 15 sec. improvement over the baseline 
H 2 /O 2 system. Further studies, both theoretical and experimental synthesis, should be carried out 
on U 3 H. 

Our previous argument concerning the stability of L 13 H would predict that the UH 3 
molecule should be unstable since back charge transfer into the H 3 ring would create an instabil¬ 
ity. In Table 4, we illustrate that this indeed is the case for UH 3 , where two separate basis sets 
both yielded a structure with imaginary frequency components, i.e. a saddle geometry. Surpris¬ 
ingly, the Li 4 structure also has an imaginary frequency for the e-mode (triangle distortion), even 
with a diffuse (6-31+ G*) basis set. The optimized geometry for Li 4 in C 3 V yields bond lengths 
which suggest that a lower order symmetry is the true ground state. A search through the 
Carnegie-Mellon University Quantum Chemistry Archive verifies this idea, where the rhombus 
structure (D 2 h) of Li 4 is found to be lower in energy. This D 2 h structure is stable, relative to 
U 3 + Li, by approximately 3.4 eV (Reference 96). Li 4 is thus probably less energetic as a fuel 
than Li(s) or LiH(s), and further studies of this molecule, considering the detrimental effect of 
molecular weight, are probably not of primary importance. 

3,2 Halogen Azides and Azide-Like Structures 

Calculations of hydroazoic acid (HN 3 ), fluorine azide (FN 3 ), fluorine isocyanate (FNCO). 
asymmetric dinitrogen dioxide (a-N 202 ) and difluoraminoborane (FNBF) were carried out at 
several levels of theory to determine the most stable geometries of these azides and azide-like 
compounds. The results are summarized in Table 5 which indicates that a stable C s structure 
exists for each of these compounds. A vibrational frequency analysis has also been carried out on 
each of these compounds. The results indicate a stable structure in 1 A' symmetry. 

A comparison with experimental structures is shown in Table 5 for HN 3 and FN 3 . It is 
clear that very good agreement has already been achieved at the MP2 level of theory. Further 
studies are in progress at higher levels of theory, and with basis sets larger than the 6-31G* sets 
used in these initial studies. 
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Chamber pressure = 800 psia 
A/A* = 20:1 
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Figure 5. Li 3 H/H 2 /F 2 Vacuum Specific Impulse 
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Table 4 . Ab Initio Structure and Calculated Energy 
for L14 and L1H3 in C3 V Symmetry 


Theory h r Energy Frequencies 






e 

a.\ 

e 

a, 

SCF/3-21G 

3.1300 

2.8191 

-29.529616 

218 i 

193 

285 

356 

SCF/3-21 + G 

3.1304 

2.8205 

-29.530374 

214/ 

192 

284 

354 

SCF/6-31G* 

3.1439 

2.8193 

-29.722764 

239/ 

199 

281 

351 

SCF/6-31 + G* 

3.1526 

2.8257 

-29.724402 

225/ 

195 

280 

349 


L1H3 


SCF/3-21G** 

1.5967 

1.2216 

-8.832136 

Unstable (2) 

SCF/6-31G** 

1.5954 

1.2294 

-8.887511 

Unstable (2) 


Energies (hartrees), distances (A), frequencies (cm -1 ) 
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Compound 


Theory 


Table 5. Optimized Geometries of Azides 



a (deg) 


0 (deg) 


hn 3 

SCF 

1.U055 

1.2381 

1.0987 

108.181 

173.815 

MP2 

1.0182 

1.2502 

1.1583 

109.946 

171.209 


EXP.l 

1.012 

1.240 

1.134 

112.65 

(180.0) 

fn 3 

SCF 

1.3820 

1.2536 

1.0995 

104.315 

174.108 

MP2 

1.4309 

1.2799 

1.1521 

103.765 

171.803 


EXP. 2 

1.444 

1.253 

1.132 

103.8 

170.9 

FNCO 

SCF 

1.3737 

1.2387 

1.1354 

109.846 

173.235 


MP2 

1.4185 

1.2622 

1.1765 

110.717 

168.914 

a-N 202 

SCF 

1.7574 

1.2024 

1.0844 

103.966 

179.506 


[MP2 

1.5305 

1.2272 

1.1548 

103.591 

179.488] 


CISD 

1.5817 

1.2240 

1.1072 

102.867 

179.344 

FNBF 

SCF 

1.2993 

1.2068 

1.2871 

180.000 

180.000 


MP2 

1.3394 

1.2491 

1.3069 

154.650 

169.291 


( ) Assumed 

[ ] V'o not stable relative to rotation to complex form 

1 Reference 118 

2 Reference 98 
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3.2.1 a-N 2 0 2 


The most significant result of our studies of azide-like structures is the prediction of a 
stable, asymmetric C s structure of dinitrogen dioxide with a characteristic IR frequency of 1206 
cm ' 1 , corresponding to the N-0 stretch. This frequency has previously been reported by Milligan 
and Jacox (Reference 97) in argon matrix studies of irradiated N 2 0 in the presence of alkali 
atoms and assigned to a possible N 2 0 2 anion. It is not characteristic of other known nitrogen 
oxides. * he calculated geometry for a-N 2 0 2 is very similar to the structures of the isoelectronic 
species FN 3 and FNCO, which have recently been characterized both theoretically and experi¬ 
mentally (References 98 and 99). 

The reaction of 0[ 1 D] atoms with C0 2 is believed to yield the cyclic C 2v CO 3 structure 
(dioxirane), which has been found by Pople et al. to be more stable than the C s structure 
(Reference 100). Therefore, calculations were also performed on the analogous C 2v structure of 
N 2 0 2 . Although it is a saddle point at the HF/3-21G level of theory, our HF/6-31G* calculations 
indicate that there is a vibrationally stable structure of N 2 0 2 about 160 kj/mol above the C, 
structure. An optimized HF/6-31G* transition state for the to C s rearrangement was found 
only 2 kj/mol above the structure. Preliminary studies at the CISD level suggest that this 
barrier becomes even smaller when electron correlation effects are included. Therefore, it seems 
unlikely that the N 2 0 2 structure will be formed in the Of'D] + N 2 0 reaction. 

There are also low-lying triplet states of N 2 0 2 , arising from 0[ 3 PJ + N 2 0 and 
0 2 [ 3 Ig] + N 2 , and it is necessary to know their location to assess the stability of the singlet. A 
CISD/6-31G* energy calculation on the lowest 3 A” state, performed at the optimized geometry 
of the singlet, shows that the triplet lies 67 kj/mol higher than the singlet. Since asymptotically the 
triplet surface lies lower than the singlet, further calculations are underway to characterize the 
crossing region and its effect on the stability of singlet a-N 2 0 2 . 

A series of Born cycle calculations were also carried out to determine the heat of formation 
of a-N 2 0 2 . The HF/6-31G result is AHf = + 430.5 kjoule/mol. Additional studies were made 
at the HF/ 6 -MG*, MP2/6-31G’, MP3/6-31G* and CISD/6-31G* levels of theory. The CISD 
calculations were corrected for size-consistency errors. The results are summarized as: 

Level of Theory AHf'(a-N 2 Q 2 ), kjoule/mol 

HF/6-31G* +559.3 
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MP2/6-31G* 

+ 493.8 

MP3/6-31G* 

+ 507.6 

CISD(SC)/6-31G* 

+ 504.5 


Excluding the HF result, which is predictably too high, we find an average value, 
AH(’ = + 504 ± 10 kjoule/mol, for the reaction: 


0 2 [X 3 Ig] + N 2 [X’X g * ] - a-NzO^A'] (76) 

This value for the heat of formation is for the gas phase species and must be reduced by the heat 
of sublimation (unknown). The energetics of the chemical synthesis of a-N 2 0 2 are shown in Fig¬ 
ure 7. Our calculated results for the equilibrium structure of N 2 0 2 are given in Table 6. 

The transition state of a-N 2 0 2 corresponds to the activation energy required to decompose 
a-N 2 0 2 into 0 2 [a‘Ag + b‘X g + ] + N 2 [X 1 S g + J. The geometry corresponding to the transition state 
is: Roo = 13980A, R N o = 1.4585A, R NN = 1.0745A, with a calculated energy of Exs 
(HF/6-31G*) = - 258.3244471 hartrees. This corresponds to an activation energy of +63.5 
kjoule/mol which is typical of all the other known azide structures. The a-N 2 0 2 dissociation 
pathway is illustrated in Figure 8. 

In addition to these thermodynamic studies, an SCF stability analysis of asymmetric 
dinitrogen dioxide (a-N 2 0 2 ) was carried out. An analysis of the stability of the RHF/6-31G* 
wavefunction for the equilibrium structure of a-N 2 0 2 has shown it to be unstable with respect to 
a broken-symmetry UHF state. Using symmetric and antisymmetric combinations of the RHF 
HOMO and LUMO for the UHF guess, a broken-symmetry UHF solution was found at about 30 
millihartrees lower energy than the RHF solution. “This state showed severe spin contamination 
(< S 2 > = 0.851). A natural orbital decomposition of the UHF wavefunction revealed two 
orbitals with occupation numbers of 1.395 and 0605, which were used as the pair orbitals in a 
GVB(l) calculation. This reorganization of the wavefunction gives an energy about 10 mil¬ 
lihartrees below the RHF and retains pure singlet character. Geometry optimizations with the 
GVB(l)/6-31G* wavefunction have led to dissociation of the 0-0 bond, suggesting that further 
investigation would be required to determine the correct zero-order wavefunction for a-N 2 0 2 . A 
limited MCSCF/6-31G* calculation for a-N 2 0 2 , which correlated all three bond pairs (O-O. 
O-N, N-N) does yield a stable structure with geometry similar to that found in the CISD/6-31G* 
calculation in Table 6. Further studies of this system are clearly indicated. 
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Available energy, kjoule/mol 




Figure 7. Energetics of the Chemical Synthesis 
of Dinitrogen Dioxide (a-N 2 0 2 ) 
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Table 6 . N2O2 Theoretical Predictions 
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Bond lengths are in angstroms, bond angles are in degrees, energies are in hartrees. and vibrational frequencies are in wavenumbers. 
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Dissociation Pathway 











Finally, in Figure 9. we show the calculated vacuum specific impulse obtainable with a 
H 2 /a-N 2 0 2 propellant combination. The peak I sp is ~ 455 sec., based on our present thermody¬ 
namic data. Metal loading would improve this performance and various hybrid combinations with 
a-N 2 0 2 are presently being studied. As a direct replacement for ammonium perchlorate (AP) in 
a typical composite propellant, we find a specific impulse improvement of 1 sec/% replacement of 
AP by a-N 2 0 2 . Using typical binder loading, a 10 second improvement should be realizable with 
current mixes. 

The molecule a-N 2 0 2 appears to be a very attractive candidate as an advanced oxidizer 
and a patent disclosure has been initiated at UTRC for both the chemical formula and several 
proposed synthesis routes for this energetic species. Attempts to synthesize a-N 2 0 2 by irradia¬ 
tion of N 2 0 in the liquid phase have been reported by Pritt (Reference 101). We have developed 
a route for synthesizing a-N 2 0 2 using the reaction of oxygen atoms, prepared in their lowest 
[*D] electronic state, with N 2 0 . In the following reactions, the molecules N 2 0 2 , N 2 and O 3 
are in the lowest or ground electronic state, O atoms are in the first excited ['DJ electronic state 
and the 0 2 molecules are in the first excited [a l A g ] electronic state. 

Oxygen atoms may be prepared in their lowest ['DJ electronic state by photolysis of N 2 0 : 


N 2 0[X 1 2 + J + hv- N 2 (X'l g + ] + 0['D], 


(77) 


or by the photolysis of O 3 : 


OjJX’A’) + hv — 0 2 [a' A g j + Of'D], (78) 

The reaction is carried out in the presence of a nonreactive collisional species, e.g. He, Ne or 
N 2 , according to the reaction: 

0[’D] + N 2 0(X’I + ] + M(E,) -* O - O - N - N[>A'J + M(E 2 ), (79) 

where M = nonreactive collisional species, and E 2 - E\ = stabilization energy of O-O-N-N in 
the lowest [ 1 A'] state. The function of the nonreactive collisional species is to kineticaily carry off 
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Chamber pressure = 1000 psia 
Exit pressure = 14.7 psia 
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Figure 9. H 2 /e-N 2 0 2 Vacuum Specific Impulse 






excess energy and allow the formation of a stable product. Reaction of the [*D] oxygen atoms 
with the N 2 O molecule would preferably be carried out in the the gas phase, however, liquid 
phase or in situu solid state formation are also possible. 

Alternatively, a - N 2 O 2 could be formed by an exchange reaction of a type: 


M'O + N 2 0[X ! X + ] -O-O-N-N^A'] + M\ 


(.SO) 


where M'O is any species in which a [ a D] oxygen atom is chemically bound. Examples of suitable 
M'O species include: HC10,HF0. ONTF, ONC1 and group fIA metal oxides where M' = Mg, Ca, 
Sn, or Ba. 

Finally, the a - N 2 O 2 product is condensed to a storable form by cryogenically cooling the 
products of the gas phase reaction or by forming a - N 2 O 2 in situa on a surface held at cryogenic 
temperatures. The condensed product collected and stored at reduced, preferably cryogenic, 
temperatures. 

Recently, Ffinchen (UTRC) has attempted synthesis by matrix isolation of the products of 
laser dissociated N 2 O. A unique IR signature has been identified in the 1220-1230 cm ' 1 region, 
very close to our predicted N-O (a ) stretch frequency. Further low temperature studies of 
a-N 2 C >2 are currently in progress to characterize this observed low temperature spectra. 
Additional theoretical studies should be carried out including characterization of the potential 
energy surface of a-N 202 in the N 2 0( 1 X'''] + 0['D] reactant region, studies of singlet-triplet 
surface crossings and studies of the stability of a-N 202 against further reaction with oxygen 
atoms or molecules. The structure and energy of possible dimers of a-N 2 02 should also be 
examined. 


3 . 2.2 FN 3 

The FN 3 system is a FIEDM material but it is also a possible chemical source for 
NF[a'A]. NFJa'A] has the potential application of pumping the BiF short wavelength chemical 
laser. Calculations of the potential energy surfaces for FN3, both ground and excited states, have 
been initiated to compliment the experimental studies of this system currently being undertaken 
at Rockwell Science Center. The energetics and pathways for decomposition, and the prediction of 
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emission/absorption wavelengths are being studied both theoretically and experimentally. These 
studies were initiated by the needs of the Air Force Weapons Laboratory (AFWL) in 
Albuquerque to model the chemistry of FN 3 decomposition. Accordingly, the scope of this 
contract was expanded to include an in-depth study of the FN 3 system. 

A detailed analysis of FN 3 has been carried out at geometries relevant to: a) the equilib¬ 
rium ground 1 A' state, b) the transition state for decomposition to NF[a'A + b'l + | 

+ N 2 [ 1 I g + ], and c) the region of singlet-triplet crossing. The details of these three studies are 
outlined below: 

a) Equilibrium Calculations 

Ab irutio theoretical calculations of the equilibrium structure and properties of fluorine azide 
(FN 3 ) have been performed using several standard basis sets and theoretical methods. SCF and 
MP2 calculations have been performed using the GAUSSIAN 86 and CADPAC programs. 6 -in -6 
CAS-MCSCF calculations have been performed using the GAMESS program. Gradient-opti¬ 
mized structures are reported in Table 7, and harmonic vibrational frequencies in Table 8 . Also 
given for comparison are the recent experimental results of Christen, Mack, Schatte, and Willner 
(Reference 98). Good agreement with the experimental structure is found, with the correlated 
methods tending to correct the short SCF bond lengths. The harmonic vibrational frequencies are 
well described by MP2, except for the highest frequency mode, the terminal N-N stretch, which is 
predicted to be over 300 cm -1 too large. Systematic studies have shown MP2 vibrational frequen¬ 
cies are usually within 5% or so for small molecules (Reference 41). The 6 -in -6 MCSCF calcula¬ 
tion reduces this error by about 100 cm' 1 , suggesting that the use of a single reference 
configuration is inadequate. 

Physically, the short terminal N-N bond length and the observed dipole moment suggests the 
electronic structure of FN 3 may be described as 

. F -6 +6 

. N - N = N : covalent form 

The MCSCF orbitals, Cl coefficients and Mulliken populations suggest that admixtures of the 
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Table 7 . FN3 Equilibirum Geometry 


N 

2 



Level 


Alpha Beta NF NN, N,N 2 Energy Dipole 

Moment 


RHF/6-31G* 

104.315 

174.108 

1.3820 

1.2536 

1.0995 

-262.602488 

1.60 

RHF/DZP 

104.345 

174.592 

1.3760 

1.2638 

1.1024 

-262.661912 

1.79 

MP2/6-31G* 

103.765 

171.803 

1.4309 

1.2799 

1.1521 

-263.282402 

1.68 

MP2/DZP 

103.836 

172.133 

1.4286 

1.2928 

1.1593 

-263.385697 

1.06 

MP2/TZ2P 

104.047 

171.324 

1.4268 

1.2695 

1.1382 

-263.601495 


MC6/6-31G* 

101.674 

174.285 

1.4630 

1.2751 

1.1154 

-262.724363 

1.99 

exp 

103.8 

170.9 

1.444 

1.253 

1.132 


1.3 

Bond lengths are in angstoms, bond angles 

are in degrees, energies are in hartrees and dipole moments 

are in Debyes. 









Table 8. 

FN3 Equilibrium Vibrational Frequencies 



Level 

a' 





a" 

ZPE 


RHF/6-31G* 

281 

757 

1045 

1225 

2386 

606 

3150 

RHF/DZP 

284 

751 

1027 

1207 

2350 

5% 


MP2/6-31G* 

247 

690 

948 

1153 

2406 

494 

2969 

MP2/DZP 

247 

677 

930 

1131 

2383 

471 


MC6/6-31G* 

241 

649 

833 

1021 

2281 

520 

2773 

exp 

241 

658 

876 

1090 

2044 

504 

2706 


Vibrational frequencies and zero point energies (ZPE) are in wavenumbers. 
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structures 


F ’ +6 -6 

■■ \ 

N = N = N ionic form - 1 

+ <5 +6 -26 
F -NsN-N: 

are significant in the MCSCF wavefunction. These would 
slightly and therefore lower its frequency. 

b) Transition State Calculations 

Calculations have been performed to locate the transition state for the decomposition of 
FN 3 into N 2 and NF. Using the GAUSSIAN 86 programs, transition state gradient optimiza¬ 
tions were performed at the RHF/6-31G* and MP2/6-31G* levels of theory. To ensure that 
stationary points of the correct curvature have been found, harmonic vibrational frequencies have 
been computed for the transition structures. The results are given in Tables 9 and 10 below. At the 
MP2/6-31G* optimized transition state, the RHF wavefunction is unstable with respect to a 
broken-symmetry UHF solution, and therefore the MP2 results may not be meaningful (for 
example, the extended central N-N bond length). Multireference calculations are required to 
reliably describe the effects of electron correlation on the transition state. Therefore, the transition 
state was also located using the 6-in-4> CAS MCSCF wavefunction. Results for the 
FN 3 -* N 2 + NF barrier height are given in Table 11. MP3 and coupled cluster calculations were 
performed at the SCF geometries to assess the effects of a more complete treatment of electron 
correlation on the barrier height. These results, included in Table 11, suggest a barrier height of 
about 0.5 - 0.7 eV, including correction for vibrational zero-point motion. 

c) Singlet-Triplet Splitting 

The location of the FN 3 1 A' -3A” crossing determines the singlet-triplet distribution of 
the FN 3 -* N 2 + NF products. Experimental evidence (Reference 102) suggests that the crossing 
is on the product side of the barrier. Ab initio SCF calculations performed to predict this crossing 
point from theory are given in Table 12. These are UHF triplet calculations done at the 1 A’ 
RHF/6-31G* optimized equilibrium and transition states. The UHF spin contamination was low. 
indicating that UHF is an acceptable level of theory for an initial approach to this problem. In 
order to locate the singlet-triplet crossing, the dissociative reaction path (intrinsic reaction 

coordinate (Reference 92)) from the RHF/6-31G* 1 A' transition state was computed using the 
GAMESS code. At selected points on the reaction path, the singlet-triplet gap was calculated. 


ionic form - 2 
(distributed) 

tend to weaken the terminal N-N bond 
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Table 9. FN 3 -» N 2 + NF Transition State 


Level 

Alpha 

Beta 

NF 

NN, 

nn 2 

Energy 

RHF/6-31G* 

169.83 

79.89 

1.3618 

1.5890 

1.0774 

-262.585108 

MP2/6-31G* 

166.87 

71.08 

1.3590 

1.8165 

1.1321 

-263.233414 

MC6/6-31G* 

171.90 

99.60 

1.3735 

1.5711 

1.1001 

-262.703604 


Bond lengths are in angstroms, bond angles are in degrees and energies are in hartrees. 


Table 10. FN 3 -► N 2 + NF TS Vibrational Frequencies 


Level 

a' 





a" 

ZPE 

RHF/6-31G* 

561 i 

186 

610 

1114 

2690 

295 

2448 

MP2/6-31G* 

410 i 

104 

427 

1077 

2113 

144 

1933 

MC6/6-31G* 

417 i 

180 

610 

1085 

2463 

290 

2314 


Vibrational frequencies and zero point energies (ZPE) are in wavenumbers. 
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Table 11. FN 3 -* N 2 + NF Barrier Height 


Level Classical (eV) Quantum (eV) 

(includes correction for 
zero point energy) 


RHF/6-31GV/RHF/6-31G* 

MP2/6-31GV/MP2/6-31G* 

MC6/6-31G *//MC6/6-31G * 

MP2/6-31GV/RHF/6-31G* 

MP3/6-31GV/RHF/6-31G* 

CCD/6-31GV/RHF6-31G* 

ST4CCD/6-31G*//RHF/6-31G* 


0.473 

0.386 

1.333 

1.205 

0.565 

0.508 

1.079 

0.999 

0.878 

0.790 

0.790 

0.710 

0.826 

0.746 


Table 12. FN 3 3A " UHF/6-31G* Energy 


Hartrees eV (Relative to ’A ) 


(w 1 A' eq -262.531687 1.927 

(w 'A' ts -262.545538 1.077 


68 





These results are shown in Table 13 and Figure 10. As SCF methods will underestimate the singlet 
energy relative to the triplet, these results should give a lower bound to the actual splitting. 
M4>ller-Plesset and coupled cluster calculations were performed at the SCF saddle point geometry 
to estimate the effect of electron correlation. These results, shown in Table 14, indicate a singlet- 
triplet splitting of about 1.8-1.9 eV 

Our results contradict those reported by Brener, et al. (Reference 103) who examined the 
singlet-triplet splitting using the MP2 optimized transition state. This state has an extended 
central N-N bond length, due to the UKF instability mentioned above, and therefore has the 
crossing on the wrong side of the transition state. The analogous singlet-triplet crossing in HN 3 
has recently been determined by Alexander and Dadigian (Reference 104) using MCSCF-C1 
methods. They find the crossing to occur at a N-N separation of 1.9 A, which is gntly larger 

th, n the 1.68 A given by our calculation. On inclusion of electron correlation effects, our 1.68 A 
separation would increase, due to the larger singlet correlation energy. We may estimate this effect 
by displacing the curves in Figure 9 to the 1.88 eV gap given by our CCD catenations. When this 

is done, the singlet-triplet crossing distance is about 1.78 A . close to the N-N distance predicted 
by Alexander and Dadigian for HN 3 . 

Additional studies of FN 3 are needed at a correlated level of theory to define .. sid 

reaction path leading to F[ 2 P] + N 3 [ 2 n] reactants. The nature of any barrier in this direction and 
an analysis of the overall thermodynamics of FN 3 formation should be carried out. 

3.2.3 CIN 3> BrN 3 

Our studies to date of FN 3 and a-tyC^ to date have illustrated that ad initio calculations 
at a correlated level of theory are extremely useful in predicting molecular structure, correlating 
experimental properties and in identifying the spectra of such species. FN 3 is now a well- 
characterized molecule and is synthesized routinely in several laboratories. This favorable 
situation, however, does not carry over to the other halogen azides where only limited spectro¬ 
scopic studies have been performed. 

Decomposition studies of halogen azides have been carried out experimentally by several 
laboratories. Gipstein and Haller (Reference 105) studied the room temperature decomposition of 
FN 3 and observed N 2 F 2 and N 2 as products. Coombe, et al. (Reference 106) examined the 193 
and 249 rm photodissociation of CIN 3 , which yields many excited singlet and triplet products. 

The thermal decomposition of CIN 3 has been studied by Richardson and Setser (Reference 107), 
who infer that NCI [ 3 Z _ ] is formed through a crossing of singlet anu triplet potential surfaces. 
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Table 13. FN 3 SCF IRC (Reaction Path) 


IRC Point 

s 

N-N 

EC A') 

E( 3 A") 

0 

0 

1.589 

-262.585108 

-262.545538 

2 

0.430 

1.663 

-262.586071 

-262.578265 

5 

0.787 

1.724 

-262.587979 

-262.601709 

10 

1.442 

1.831 

-262.592846 

-262.635125 

20 

2.778 

2.047 

-262.603338 

-262.678595 


Bond lengths are in angstroms and energies are in hartrees. 


Table 14. FN 3 1 A' - 3A" Higher Order Energies at SCF Geometry 


Level 


'A' 

RHFts 


3 A" 

@ RHF Is 


Energy 

Difference (eV) 

3 A" - 'A' 


SCF/6-31G* 

-262.585108 

-262.545516 

1.077 

PUHF 


-262.551575 

0.912 

MP2/6-31G* 

-263.221058 

-263.148967 

1.962 

PUMP2 


-263.153153 

1.848 

MP3/6-31G* 

-263.218597 

-263.149187 

1.889 

CCD/6-31G* 

-263.225497 

-263.156454 

1.878 


Energies are in hartrees. 
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TOTAL ELECTRONIC ENERGY (HARTREES) 









BrN 3 decomposition has been studied by Coombe (Reference 108) who suggests that the N-N 
bond in BrN 3 is weak, relative to the other hydrogen azides. 

Benard and associates (Reference 109) (Rockwell Science Center) have studied the 193 nm 
photolysis of FN 3 and observed both NFfb 1 ^) and NFfa'AJ product states. Their work 

suggests that either excited F^f'A'j surfaces are involved or that a series of singlet/triplet 
crossings must occur. More recently. Benard. et al. have carried out a unique thermal dissociation 
experiment with FN 3 which indicates that, with thermal energy transfer, only the NFfa'A) state 
is formed. Their adiabatic decomposition mechanism correlates with the location and barrier 
height predicted by our theoretical studies of this system. 

It is clear from the above that the stability and decomposition mechanisms of the halogen 
azides need to be systematically studied. The key points that need to be addressed theoretically 
are: 


1) heat of formation of XN 3 (X = F, Cl, Br) 

2 ) zero-point energy analysis 

3) vibrational (IR) spectra of the ground 1 A' state 

4) barrier height and width of the transition state leading to NX + N 2 

5) location of low-lying singlet and triplet reaction surfaces 

6 ) analysis of singlet-triplet crossings and branching into NX product channels. 

As mentioned above, our current theoretical studies of advanced oxidizer types under this 
research program have included a preliminary study of the FN 3 molecule. In collaboration with 
Dr. David Benard’s experimental group at Rockwell Science Center, we have made significant 
progress in characterizing the thermochemistry and stability of FN 3 . Our calculated barrier 
height (—0.6 eV) for decomposition of FN 3 ['A') into NFfa'A] + N 2 [X 1 I g + ) and our predicted 
heat of formation ( + 565 kJ/mol) have both been confirmed in Benard’s studies On the basis of 
our past experience, we believe that collaborative efforts are especially fruitful in a systematic 
characterization of the chemistry of the halogen azide series. Calculations of reaction barrier 
heights could be used to guide experimental dissociation studies. Experimental frequencies could 
also be correlated against the prediction of electronic structure calculations. 

3.3 Cyclic Boron Structures 

In addition to our studies of Li m H n structures, a preliminary analysis of several beryllium 
and boron compounds was undertaken. There are several beryllium-containing structures that 
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appear promising. Linear BeH 2 is well characterized both theoretically and experimentally 
(Reference 110). However, BeHi in C s symmetry, which is essentially a BeH + + H“ ionic 
complex, has not previously been considered and was a candidate for theoretical study. Nicolaides 
and Zdetsis (Reference 23) have examined the corresponding HeH 2 (C s ) molecule but no stable 
ground state is indicated other than He + H 2 . The molecular structure, LiBeH 2 (C 2 V ), was also 
suggested since the Be atom is very electron deficient in BeH 2 and a Lewis adduct may be 
possible with an alkali. Unfortunately, no stable structure could be found for either of these 
species, other than the well-characterized D»h linear structure of BeH 2 . The performance of 
B c H 2 (D* h ) as a fuel has previously been analyzed ( Reference 111) but its potential is limited by 
environmental considerations based on the toxicity of BeO. 

Beyond beryllium, the potential for forming high energy fuels, which offer a significant 
increase in I sp , diminishes rapidly owing to the molecular weight considerations discussed above 

and illustrated in Figure 1. However, there are several boron compounds that might be promising 
candidates, provided they represent local regions of stability and offer a possible route of 
synthesis. Diborane (B 2 H 6 ) has been extensively studied and evaluated as a rocket fuel, the 
maximum specific impulse ( P c = 1000 psia, expansion to sea level) available with O 2 as an 
oxidizer is 342 seconds, a value below that offered by H 2 /O 2 (389 sec) and Be/Hi/C^ (459 sec). 
The detrimental effects of molecular weight are apparent with pentaborane (B 5 H y ) and decabo- 
rane (Bi 0 Hi 4 ) , which yield 319 and 282 seconds, respectively. This led us to suggest that BH 3 
adducts may offer higher energy, and hopefully may represent isomeric structures with some 
measure of stability. Some other possibilities here, which have been examined theoretically by 
Lammertsma’s group, are the HjB-Li adducts and possibly a H 3 B-Be structure, similar to 
the proposed BeH 2 (C,) molecule. 

There are several other possibilities if we consider B 2 or B 3 structures. A few compounds 
which are suggested from similar hydrocarbon structures are shown in Figure 11. Cvclotriborane 
(B 3 H 3 ) would be the saturated equivalent of cycloproplyene (C 3 H 4 ) in hydrocarbon chemistry. 
Such a highly strained structure might form substitution compounds which have greater stability 
but would still be of value as potential fuels. Iminodiborane (B 2 H 2 NH), which is analogous to 
aziridine (ethylencimine), would have a molecular weight of only 38 and possibly a very high 
positive heat of formation. Its structure (or stability) was totally unknown. Another possibility is 
1,2 epoxydiborane (B 2 H 2 O), the equivalent of ethylene oxide in hydrocarbon chemistry. A 
somewhat heavier (MW = 57) triborylineimine or iminotriborane (B 3 H 3 NH). which is similar to 
azetidine (trimethyieneimine) might also be promising. Similar highly strained compounds 


73 




Cyclotriborane (B 3 H 3 ) 


H 



O 0 

R BB = 1,734 A R b -H' = 1.180 o A R b . b = 1.536 o 

R BH = 1,178A R bh = 1.174A R bb = 2.032 A 

e = - 75.174458, unstable (2) E = - 75.647600, stable 

Iminodiborane (B 2 H 2 NH) 

H 

I 


N 



R BB = 1.621 A L HBN = 136.2° 

R BH = 1,174A ^ BNH = 145.3° 

R nh = 0.989 A E= - 104.981636 

R BN = 1,423 A 


Frequencies (cm - "*) 



700 818 921 923 989 1000 1110 1268 1423 2878 2916 3965 


Figure 11. Cyclic Boron Compounds 
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containing four or more boron atoms would appear to offer little hope of yielding a specific 
impulse above 400 seconds. 

In a preliminary study of cyclic boron structure, under the present contract, a limited set of 
calculations were carried out on cyclotriborane (B 3 H 3 ) and iminodiborane (B 2 H 2 NH) at the 
SCF level of theory using a 6-31G* basis set. The C 3 V structure for B 3 H 3 is found to be unstable 
but further studies in C 2 V symmetry indicate stability. The anion B 3 H 3 is found to be stable as 
predicted by Lipscomb’s rules. Iminodiborane (B 2 H 2 NH), on the other hand, is a stable 
structure in C 2 v symmetry. The optimized geometries and calculated frequencies are shown in 
Figure 11. Further studies of B 2 H 2 NH are indicated but our calculated energetics to date 
indicate that this is not a high-energy structure. We find AH? — - 20 kcal/mol for B 2 H 2 NH , 
which can be compared with AH? = + 12.3 kcal/mol for diborane. A further negative view is 
that most iminoboranes eventually form borazin'. (B 3 N 3 H 6 ), which is not a high energy fuel. 
Further studies of other cyclic boron compounds, however, may prove useful. 

3.4 Hypervalent Structures 

A third category of compounds which may represent energetic molecular conformations are 
those resulting from orbital hybridization to a high valtnce state of the central coordinating atom 
in the molecular species. The existence of such species and the concept of hypervalency have been 
described many years ago by Pauling (Reference 112). The successful synthesis of hypervalent 
compounds of second-row elements, such as PF 5 , H 3 PO, H 3 PF 2 , IF 7 , XeFfi, has been carried out 
in several laboratories. The textbook explanation of bonding in such compounds involves 
contribution from the normally unoccupied d-orbitals on the coordinating atom to form sp*d 
(square planar), sp 3 d (trigonal bipyramidal), or sp 3 d 2 (octahedral) structures. An extension of 
this argument is that hypervalent molecules formed from first row atoms, such as NF 5 , are not 
possible since the hybridization energy is too large. 

In recent years, quantum mechanical studies have revealed that the bonding in a hypervalent 
structure such as PF 5 is much more complicated than the simple Pauling hybridization model. In 
fact, detailed calculations (Reference 113) reveal that PF 5 is stable, even if the d-orbital partici¬ 
pation is set to zero. The bonding appears to be better described by sp 2 hybridization to form 
the trigonal PF 3 ring, with the axial fluorines in a 3-center sp bond. Such a structure can be 
further stabilized by back-bonding contributions of the normally unoccupied d—orbitals. This new 
insight into the bonding of hypervalent compounds revives the possibility of forming such 
structures from first row atoms. Compounds such as NF 5 , NH 3 F 2 , OF*. CIF 5 O and CIF 5 
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represent hypervalent, high oxidation state structures which represent a new class of possible 
oxidizers for chemical propulsion systems. For example, Christe (Reference 114) has estimated 
that theoretical performance calculations for H 2 /NF 5 yield a vacuum I sp of 505 sec. We have 

estimated the optimum performance for a H 2 /F 2 /NF 5 tripropellant (P c = 800 psia, A/A * = 
20:1) to be 466 seconds, yielding no significant improvement over the H 2 /F 2 system. However, if 
NF 5 has physical properties as a storable oxidizer, it may be a promising new candidate in 
comparison with CIF 5 . 

The synthesis of such strong oxidizers is an obvious experimental challenge. Christe and 
collaborators (Reference 115) have succeeded in the synthesis of the adduct, tetrafluoronitronium 
(V) hexafluoroarsenate (NT 4 AsF,>) but all attempts to date to synthesize NF 5 have failed. The 
structure of NF 4 AsF* has been confirmed by spectroscopy and magnetic resonance measure¬ 
ments and is known to be ionic, NF 4 + AsF* . The hybridization of the N atom in NF 4 AsF a is thus 
sp 3 , forming the stable tetrahedral NF 4 + cation, rather than a true pentavalent structure. 

Further, semi-quantitative arguments by Kutzelnigg (Reference 116) suggest that NF 5 or NH 3 F 2 
structures will not exist for steric reasons; the charge distribution in the NF 3 planar ring excludes 
the formation of a short axial N-F bond. 

In an effort to more fully understand the competition between steric factors, which might 
prevent axial bond formation, and the strengths of 3-center bonds, an analysis of NF 5 structures 
was carried out at both the SCF and MP2 levels of theory using several basis sets. An ionic 
NF 4 + F~ structure was first examined in light of the stability of Christe’s NF 4 + AsFfc~ structure. 
No stable minimum could be located starting the geometry search along the NF 4 + + F* reaction 
path. This structure simply decomposes to NF 3 + F 2 . However, two stable covalently bonded 
structures were found for NF 5 at the SCF level. These structures are illustrated in Figure 12. The 
D 3 h (trigonal bipyramid) structure was found to be lower in energy at the SCF/6-31G* level of 
theory and further, MP2 calculations have indicated that the C 3v structure (inverted umbrella) 
represents a saddle point. The calculated geometry and vibrational frequencies for the D 3 h 
structure are given in Tables 15 and 16, respectively, at several levels of theory. These preliminary 
calculations clearly indicate that careful optimization of the basis set and a correlated level of 
theory will be required to positively identify the structure and vibrational modes. Introduction of 
diffuse basis functions and optimization of the d-orbital exponents have a small but important 
effect on the thermochemistry. Calculations at the MP2 level of theory clearly show a different 
vibrational pattern owing to the lengthening of the planar NF bond. The most significant 
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Table 15. Ab Initio Structure and Calculated Energy for NF 5 in Djh Symmetry 


Theory d-exponent r p r h Energy 

N F 


SCF/6-31G 

- 

- 

1.3860 

1.5926 

-550.936601 

SCF/6-31G* 

0.80 

0.80 

1.3267 

1.5312 

-551.103856 

SCF/6-31G(d) 

0.94 

0.76 

1.3260 

1.5310 

-551.104823 

SCF/6-31+ G* 

0.80 

0.80 

1.3223 

1.5473 

-551.126800 

SCF/6-311G 

- 

- 

1.3765 

1.5978 

-551.116124 

SCF/6-311G* 

0.91 

1.75 

1.3132 

1.5457 

-551.249831 

SCF/6-31 lG(d) 

0.94 

0.76 

1.3260 

1.5310 

-551.274677 

SCF/6-311G(2d) 

1.63 

1.32 1 

1.3148 

1.5429 

-551.292792 


0.54 

0.44 * 




MP2/6-31G* 

0.80 

0.80 

1.4098 

1.5472 

-552.228780 

MP2/6-31G(d) 

0.93 

1.47 

1.3962 

1.5538 

-552.302405 


Energies (hartrees), distances (A) 
r p = N-F bond length in the plane 
fh = N-F axial bond length 
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Table 16. Calculated Harmonic Frequencies for NF S as a D 3h Structure 


sO r- 

go r- oo 

n — — 


O O O 


r** i m 

si si si si 


V3V5C/5V5C/5<y5V2C/5 
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observation, however, is that NF 5 is predicted to be vibrationally stable without incorporation of 
d-orbitals. The classical model of sp 3 d bonding is thus clearly incorrect for this system. To check 
this observation, calculations were carried out on PF 5 at the SCF/6-31G level of theory, both with 
and without d-orbitals in the basis. Similar stable vibrational structures are found in either case 
so the energetics of hybridization to sp 3 d is not the underlying factor in the stability of these 
structures. 

At the SCF level, the a. 2 " stretch is the lowest frequency vibrational mode for NF 5 . 
Decomposition along this normal mode will probably lead to NF 4 + F, an endothermic step. 
However, at the MP2 level, the e' axial bend is predicted to have the lowest frequency. Decompo¬ 
sition along this normal mode would lead to NF 3 + F 2 , clearly an exothermic process. Calcula¬ 
tions for PF 5 also indicate that e' is the lowest normal vibrational mode. Detailed calculations of 
both equilibrium geometries and transition states are clearly required here. 

In an attempt to understand the dependence of the a 2 ” (S 3 ) mode of NF 5 on the diffuse¬ 
ness of the basis, several calculations of NF 5 were carried out in C 3 V symmetry, corresponding to 
a displacement along the S 3 coordinate. We find a C 3 V stationary point at the SCF/6-31G* level 
of theory which is only 0.1 eV higher than the corresponding D 3 h structure. A harmonic 
frequency analysis indicates that the C 3 V structure has all real frequencies with the a' (3) axial 
stretch mode corresponding to the weakest displacement. This observation suggests that NF 5 , at 
this level of theory, has a small barrier along this coordinate displacement between the nearly 
isoenergetic and stationary D 3 h and C 3 V structures. 

To test the stability of the C 3 V structure at a higher level of theory, MP2 calculations were 
carried out with the 6-31G* and 6-311+ G* basis sets. The starting geometry corresponded to the 
optimized C 3 V structures found at the SCF level of theory. In both cases, the optimization 
proceeded uniformly to the lower energy D 3 h structure. A comparative examination of the SCF 
charge distributions in the C 3 V and D 3 h structures indicates that the C 3 V structure is very ionic 
[NF 4 d F' d ] relative to the more covalently bonded D 31 , structure. The more diffuse basis sets 
favor this ionic distribution, resulting in the indicated stability at the SCF level. At the MP2 level 
of theory, the more ionic C 3 V structure lies higher in energy and geometry optimization always 
yields the lower energy, and more covalent, D 3 h structure. 

A final point can be made concerning the possibility of synthesis of this compound. 

Although a definitive study of the thermodynamics of formation of NF 5 has not yet been 
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completed, it is clear that the formation of NF 5 is endothermic relative to NF 3 + F 2 , whereas this 
process is exothermic in the case of PF 5 . It also appears unlikely that NF 5 can be synthesized 
starting from a salt-type structure NF 4 + R~, since the known tetrahedral compounds have 
correspondingly lower heats of formation. The key to successful synthesis may lie in the initial 
preparation of the reactants NF 3 and F 2 . By expanding the bond lengths in NF 3 and F 2 
through vibrational excitation, such that the axial and equatorial N-F bonds are more nearly 
equal, the insertion of F 2 above the NF 3 equatorial plane, followed by a Berry psuedorotation 
(Reference 117) of an equatorial F atom into an axial position below the plane may occur. This 
same insertion mechanism is believed to occur in the synthesis of PF 5 . Rapid collisional cooling 
of the vibrationally hot NF 5 may result in trapping of a stable product. 
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SECTION 4 


CONCLUDING REMARKS 

The main goal of this research investigation was to identify and characterize several high 
energy species that have potential for advanced chemical propulsion systems. Using deliverable 
specific impulse and storability as the major criteria for the evaluation oi new compounds, 
primary consideration was given to ground state molecular structures, of low molecular weight, 
which exhibit a high positive heat of formation. We can categorize the compounds that were 
studied into four (4) groups: 

1 . light element metal and metal hydride structures as advanced fuels 
(H 4 .Li 3 H, LiH 3 , Li 4 ) 

2. azide-like structures as advanced oxidizers (FN 3 , a - N 2 C> 2 , FNCO, C0 3 , HN 3 
FNBF) 

3. cyclic boron/boron hydride structures as advanced fuels ( B 3 H 3 , B 2 H 2 NH ) and 

4. hypervalent compounds as energetic advanced oxidizers ( NF 5 and for comparison, 

PF S ). 

A summary of our findings in each of these groups follows: 

1. Wc conclude that the ground states of H 4 (tetrahydrogen), LiH 3 and Li 4 are not 
promising candidates for storable fuels. None of these molecules exhibit chemical 
stability in a highly energetic conformation. The system Li 3 H, as a planar C 2 v 
structure, is stable in its ground state and may offer some promise as a component 
in a solid hybrid configuration. More work is needed, both theoretically and exper¬ 
imentally on this compound. In particular, the stability against dimer formation, 
with subsequent formation of LiH and H 2 , needs to be examined. 

2. All of the azide-structures are vibrationally stable with high positive heats of forma¬ 
tion, excepting FNCO The drawback with these materials is their high molecular 
weight and high chemical reactivity. However, both FN 3 and a-N 202 , if encapsu¬ 
lated at about the 10% level as replacement for ammonium perchlorate (AP), yield 
a significant improvement in the I sp of typical solid formulations. In particular, a 

10% replacement of AP with a-N 202 yields a 10 second improvement, about the 
same as that found by the addition of A1 metal. Experimental work on a-N 202 . 
synthesis and property evaluation, is indicated. 
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3. The cyclic boron structures B 3 H 3 and B 2 H 2 NH proved to be less energetic than 
was anticipated. B 2 H 2 NH is more stable than diborane with a heat of formation 

— - 30 kca!/mol, relative to B 2 HA . The propensity of such compounds to undergo 
condensation reactions forming borazines; and their high molecular weight, suggest 
that these cyclic structures are not prime candidates for consideration as new fuels. 

4. The only hypervalent compound studied extensively to date is NF 5 . This molecule 
has stability as a trigonal bipyramid structure, with a predicted heat of formation, 
AH°(0 K) = +25.1 kJ/mol. We have outlined a possible synthesis route, starting 
with NF 3 which has been vibrationally excited via collisional or laser pumping. This 
system has interesting possibilities as a new liquid oxidizer or in a hybrid formula¬ 
tion with F 2 or O 2 , depending on the physical properties of NF 5 . An experimental 
program to synthesize NF 5 , starting from NF 3 , is recommended, based on our ex¬ 
isting knowledge of the properties of this system. 

A final comment is that, based on these studies, we predict advances in oxidizer formula¬ 
tions should be possible, particularly in the area of new solid propellant formulations. We are 
optimistic that a new oxidizer can be developed through further theoretical and experimental 
studies which follow this initial investigation. We are not as optimistic concerning the develop¬ 
ment of a new fuel. Considering the low molecular weight and relatively weak chemical bond, it 
may prove difficult to find a better fuel than H 2 unless progress in free radical stabilization or 
entrainment of light metals in H 2 is achieved. 
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APPENDIX B-l 


ELECTRONIC STRUCTURE AND STABILITY OF SMALL 
CATION AND ANION HYDROGEN CLUSTERS 

J.A. Montgomery, Jr. and H.H. Michels 
United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

Ab initio calculations of the electronic structure of H„ and H,t clusters have been carried out 
using accurate Gaussian basis sets and with levels of theory up to fourth-order perturbation theory 
(MP4) and single and double excitation configuration interaction (CISD). The odd hydrogen cation 
addition sequence, H„ + H 2 -* H„ + 2 , appears to be thermodynamically stable for large size cluster 
formation. The even H£ cation also exhibits surprising stability in D 2 J symmetry. In contrast, the 
hydrogen anion addition sequences, H„ + H-* H „+1 and H„ + H 2 -* H„ + 2 • appear to be 
thermoneutral or unstable. 


Supported in part by APOSR under Contract F49620-85-C-0095 and in part by 
AJFRPL under Contract F04611-86-C-0071. 

Published in the Proceedings of the Cooling, Condensation and Storage of 
Hydrogen Cluster Ion Workshop, Editor - J.T Bahns, P 83, 1987. 
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APPENDIX B-2 


ON THE STABILITY OF H 4 IN C 3v SYMMETRY 


J.A. Montgomery, Jr. and H.H. Michels 
United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

The electronic structure of H 4 in Cj v symmetry has been studied by ub initio calculations. The 
ground state energy minimum found in these calculations is shown by vibrational analysis to be a 
saddle point rather than a stable bound state. The path of steepest descent from the C 3 V saddle point 
is found to lead monotonically to dissociation into two H 2 molecules. 


Supported in part by APAL under Contract F04611-86-C-0071. 
Published in Journal of Chemical Physics. 86, p. 5882, 1987. 
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APPENDIX B-3 


ON THE STRUCTURE OF THE GROUND STATE OF 


J.A. Montgomery, Jr. and H.H. Michels 
United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

The electronic structure of H* has been studied by ab initio calculations. It is found that there 
exist two stable structures, of C, and D 2 d symmetries. Post-SCF calculations suggest that the 
ordering of these states is determined by electron correlation, with the D 2 d structure lying slightly 
lower than the C s structure. It is suggested that the D 2 d structure is formed in the recent experiments 
of Kirschner and Bowers. 


Supported in part by AFAL under Contract F04611-86-C-0071. 
Published in Journal of Chemical Physics, 86, p. 5882, 1987. 
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APPENDIX B-4 


ELECTRONIC STRUCTURE AND STABILITY OF 
ENERGETIC CHEMICAL SPECIES 


H.H. Michels and J.A. Montgomery, Jr. 
United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

Theoretical quantum mechanical calculations have been carried out for several potentially high 
energy chemical species. Our research objectives are to identify those light element molecular species 
that have sufficient internal energy to be useful in advanced chemical propulsion systems. A 
parametric study of the deliverable specific impulse as a function of available reaction energy and 
propellant weight indicates that primary consideration should be given to molecular structures that 
can be formed from hydrogen through boron and that have a molecular weight of 40 or less. 

Calculations to date have been performed on three classes of light element compounds: C 3 V 
structures such as H 4 and L^H, azide structures such as N 3 F and a-N 2 0 2 and cyclic boron 
structures such as B 3 H 3 and B 2 H 2 NH. We find that the ground state potential energy minimum 
found for H 4 in C 3 ,, symmetry corresponds to a saddle region rather than a stable bound state and 
that distortion via vibrational or rotational modes leads monotonically to dissociation into two H 2 
molecules. In contrast, L^H appears to be chemically stable in C 3 V symmetry. Our studies of 
asymmetric dinitrogen dioxide (a-N 2 0 2 ) indicate a stable azide-like structure in C s symmetry for the 
ground 1 A' state, with an indicated heat of formation of +430 kjoule/mol. This molecule is a very 
attractive candidate as an advanced oxidizer. Our preliminary calculations for boron compounds 
indicate that B 3 H 3 is unstable in D 3 h symmetry but that iminodiborane (B 2 H 2 NH) is stable in 
symmetry. The electronic structure, vibrational analysis and thermodynamic stability of these 
compounds will be described. 


Supported in part by AFAL under Contract F04611-86-C-0071. 

Published in Proceedings of the High Energy Density Matter Conference, p. 219-238. 1987. 
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APPENDIX B-5 


THE ELECTRONIC STRUCTURE AND STABILITY OF 
ASYMMETRIC DINITROGEN DIOXIDE (a-N 2 0 2 ) 

H.H Michels and J.A. Montgomery, Jr. 

United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

The electronic structure of asymmetric dinitrogen dioxide (a-N 2 0 2 ) has been studied by ub initio 
calculations. A stable structure in C 5 symmetry, of similar geometry to FN 3 , is predicted for the 
lowest 'A' state. This structure is characterized by an IR frequency of 1206 cm' 1 for the N-O stretch 
and exhibits a significant activation barrier (65 kJ/mol) for decomposition into N 2 + singlet 0 2 . This 
asymmetric form of N 2 0 2 has a heat of formation of approximately +500 kJ/mol and is potentially 
explosive. 


Supported in pan by AFAL under Contract F04611-86-C-0071. 
Published in Journal of Chemical Physics, 88, p. 7248, 1988. 
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APPENDIX B-6 


THEORETICAL INVESTIGATION OF ENERGY STORAGE 
IN ATOMIC AND MOLECULAR SYSTEMS 

H.H. Michels and J.A. Montgomery, Jr. 

United Technologies Research Center 
East Hartford, Connecticut 06108 

ABSTRACT 

Theoretical electronic structure calculations are being carried out for several high energy species 
that are attractive candidates for advanced chemical propulsion systems. Using deliverable specific 
impulse and storability as the major criteria for the goodness of new oxidizers or fuels, primary 
consideration is being given to ground state molecular structures, of low molecular weight, which 
exhibit a high positive heat of formation. Calculations to date have been carried out on: 1 ) light element 
C 3 V structures H 4 , L^H, UH 3 , Li 4 ); 2) azide-like structures (FN 3 , a-N 202 , FNCO, CO 3 , 
FNBF); 3) cyclic boron structures (B 3 H 3 , B 2 H 2 NH); 4) hypervalent structures (NF 5 ). 

We find that the ground states of H 4 (tetrahydrogen), LiH 3 and Li 4 , as C 3 V structures, are 
unstable toward dissociation into diatomic fragments. These molecules do not look promising as 
storable, high energy chemical species. L^H appears to be vibrationally stable as a C 3 V structure but 
preliminary studies indicate a low energy content. 

All of the azide-like structures are vibrationally stable with high (400-500 kJ/mol) positive heats 
of formation. The best candidate appears to be asymmetric dinitrogen dioxide, a - N 2 O 2 . which 
exhibits a heat of formation of +500 kJ/mol and can deliver 450 seconds I sp as a hybrid propellant 
with hydrogen. 

None of the cyclic boron structures examined to date are highly energetic. Their chemistry is 
further complicated by condensation reactions which cause conversion to higher molecular weight 
structures. 

The hypervalent NF 5 molecule has been studied as a D 3 ;i structure, as the ionic NF 4 F“ form 
and as a C 3 V structure. Preliminary calculations indicate stability for NF 5 both as a D 3 h and C 3 V 
structure. Detailed correlated energy calculations are in progress for NF 5 and for several azide 
structures. 


Supported in part by AFAL under Contract F04611-86-C-0071. 

Published in Proceedings of the Air Force Second High Energy Density Materials Contractors 
Conference, Editors - L. P. Davis and F. J Wodarczyk, p. 93-100, 1988. 


B-6 





APPENDIX B-7 


REPLY TO THE COMMENT ON: ON THE STABILITY OF 
H 4 IN C 3V SYMMETRY 


J.A. Montgomery, Jr. and H.H. Michels 
United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

This is a reply to a response by E. M. Evleth and E. Kassab, On the Stability of H 4 in Cj v 
Symmetry: Comments, JCP, 89, p. 3928,1988. Their response was to an earlier publication by the above 
authors (JAM and HHM) entitled: On the Stability of H 4 in C 3 V Symmetry, JCP. 86, p. 5882, 1987. 


Supported in part by AFAL under Contract F04611-86-C-0071. 
Published in the Journal of Chemical Physics, 89, p. 3929, 1988. 
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APPENDIX B -8 


THEORETICAL INVESTIGATION OF ENERGY STORAGE 
IN ATOMIC AND MOLECULAR SYSTEMS 

H.H. Michels and J.A. Montgomery, Jr. 

United Technologies Research Center 
East Hartford, Connecticut 06108 

ABSTRACT 

Theoretical electronic structure calculations are being carried out for several high energy species 
that are attractive candidates for advanced chemical propulsion systems. Using deliverable specific 
impulse and storability as the major criteria for the evaluation of new oxidizers or fuels, primary 
consideration is being given to ground state molecular structures, of low molecular weight, which 
exhibit a high positive heat of formation. Calculations to date have been carried out on: 1 ) light element 
C3V and C 2 v structures (H4, LfyH, LiH 3 , Li.t): 2 ) azide-like structures (FN3, a-N 202 . FNCO, 
CO 3 , HN 3 , FNBF); 3) cyclic boron structures (B 3 H 3 , B 2 H 2 NH); and 4) hypervalent structures 

nf 5 , pf 5 ). 

We have studied the L^H system in both C 3 V (trigonal pyramid) and C 2 v (planar kite-like) 
symmetries. Both structures are vibrationally stable, with the C 2 v symmetry lower in energy at the 
MP2/6-311G** level of theory. Using either O 2 or F 2 as an oxidizer. Li 3 H appears to give about 10% 
improvement in I sp over LiH. 

The energies and pathways for decomposition of FN 3 have been studied. Our calculations, 
carried out to compliment the experimental studies of FN 3 at Rockwell Science Center, indicate that 
NF[a‘AJ should be efficiently produced by the thermal decomposition of FN 3 . The structure of this 
molecule, its predicted vibrational spectrum and the role of low-lying triplet states have been 
examined. A parallel study of the isoelectronic FNCO molecule is in progress. 

Further studies of the hypervalent NF 5 molecule indicate vibrational stability only in D 3 h 
symmetry. A comparison of the bonding in NF 5 and PF 5 is currently in progress. 

A preliminary experiment has shown evidence for formation of a - N 2 O 2 from the reaction of 
O^Dj with N 2 O in a neon matrix on a window cooled to -10K. A complete vibrational 
characterization, to positively identify that a - N 2 O 2 has been formed, is in progress. 

Supported in part by AFAL under Contract F04611-86-C-0071. 

Published in Proceedings of the Air Force Third High Energy Density Matter 
(HEDM) Contractors Conference, Eds. — T. G. Wiley and R. A. van Opijnen, 
p. 61-72, 1989. 
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APPENDIX B-9 


THE STRUCTURE AND BONDING OF Li 3 H 
ION-PAIR STATES 


J.A. Montgomery, Jr. and H.H. Michels 
United Technologies Research Center 
East Hartford, Connecticut 06108 

O. F. Guner and K. Lammertsma 
Department of Chemistry 
University of Alabama at Birmingham 
Birmingham, Alabama 35294 


ABSTRACT 

Ab initio calculations on L 13 H ion-pair states and a topological analysis of the charge density of 
the resulting optimized structures are reported. The global Li 3 H minimum is a planar C 3v structure 
with H'-Li 3 + ion-pair character. The pyramidal C 3v structure, which is 21.5 kcal/mol higher in 
energy (MP2/6-311 + + G * *), is found to have ionic H-Li bonds, bu; no Li—Li bonds. Although 
vibrational analysis indicates this structure to be a true minimum on the potential energy surface, 
inclusion of the vibrational zero-point energy makes the thermodynamic stability of this structure 
questionable. 


Supported in part by AFAL under Contract F04611-86-C-0071. 
To be Published in Chemical Physics Letter, 1989. 
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APPENDIX B-10 


THE ELECTRONIC STRUCTURE AND STABILITY OF NF 5 


H.H. Michels and J.A. Montgomery, Jr. 
United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

An ab initio study of the electronic structure of NF 5 and PF 5 has been carried out using 
M<$>ller-Plesset perturbation theory. Optimized geometries were calculated at the SCF and MP2 levels 
of theory using several basis sets, ranging from 6-31G to 6-311 + G*. A vibrational analysis indicates 
that NF 5 , in D 3 h symmetry, has all real frequencies, even if d orbital contributions are set to zero. 

Dissimilar N-F bond lengths (r eq = lJSA.ra, = 1.58A) are found, in contrast to PF 5 which exhibits 
nearly equal equatorial and axial bond lengths. A topological analysis of the calculated charge 
distribution in NF 5 indicates true pentavalent coordination. Our results suggest that synthesis of 
hypervalent compounds of first row-atoms may be possible. 


Supported in part by AFAL under Contract F04611-86-C-0071. 
To be published in Chemical Physics Letters, 1989. 
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APPENDIX B— 11 


THEORETICAL STUDIES OF THE VIBRATIONAL 
SPECTRUM AND THERMAL DECOMPOSITION OF 
FLUORINE AZIDE (FN 3 ) 


J.A. Montgomery, Jr. and H.H. Michels 
United Technologies Research Center 
East Hartford, Connecticut 06108 


ABSTRACT 

The vibrational spectrum and decomposition pathways of the highly energetic FN 3 molecule 
have been studied by ub initio methods. Using extended basis sets, harmonic vibrational frequencies 
and isotope shifts have been calculated from SCF and Mcfcller-Plesset perturbation theory for the 
ground 'A' state of FN 3 . We find good agreement with the experimental molecular geometry and 
observed IR frequencies at the MP2/6-31G* level of theory. The transition state for the thermal 
decomposition of FN 3 into N 2 and NF has been located using both MCSCF and perturbation theory- 
methods. We find a significant dependence of the calculated transition state barrier height on the level 
of theoretical description. At the MP2/6-31G* level, the optimized RHF wavefunction for the 
transition state is unstable with respect to a broken-symmetry UHF solution, and therefore 
multi-reference calculations are required to properly describe the effects of electron correlation. The 
results are compared with existing experimental spectroscopic data and thermal decomposition 
studies. 
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APPENDIX C 


THE STRUCTURE AND BONDING OF Li 3 H 
ION-PAIR STATES 

J. A. Montgomery, Jr. and H. H. Michels 
United Technologies Research Center 
East Hartford, CT 06108 

O. F. Guner and K. Lammertsma 
Department of Chemistry 
University of Alabama at Birmingham 
Birmingham, AL 35294 


ABSTRACT 

Ab initio calculations on Li 3 H ion-pair states and a topological analysis of the charge 
density of the resulting optimized structures are reported. The global Li 3 H minimum is a planar 
C 2 v structure with H'-Lij - ion-pair character. The pyramidal C 3v structure, which is 21.5 
kcal/mol higher in energy (MP2/6-311 + + G * *), is found to have ionic H-Li bonds, but no 
Li-Li bonds. Although vibrational analysis indicates this structure to be a true minimum on the 
potential energy surface, inclusion of the vibrational zero-point energy makes the thermodynamic 
stability of this structure questionable. 
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